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Statistical  Mechanics  and  Rheological  approaches  have  been  utilized 
to  establish  a  mathematical  model  to  represent  the  long  term  creep  of 
concrete. 

The  generalized  model  presented  predicts  behavior  of  lightweight  and 
conventional  concretes  irrespective  of  size  and  shape  variations.  A 
coefficient  of  structural  stability  has  been  introduced  to  incorporate  the 
effect  of  aging  and  of  stress  variations  in  the  analysis  of  data.  The 
use  of  the  age  and  stress  dependent  components  of  the  stability  coefficient 
permits  the  calculation  of  creep  at  any  time  and  for  any  stress  level. 


CHAPTER  I 


INTRODUCTION 


The  mechanical  properties  of  concrete,  and  perhaps  all  matter,  may 
be  treated  from  three  distinct  viewpoints;  each  viewpoint  attenq>ts  to 
satisfy  an  engineering  need  for  predictability  concerning  the  mechanical 
behavior  of  the  material  in  specific  structural  applicaitions;  each 
viewpoint  used  to  satisfy  this  need  takes  a  basically  separate  route; 
each  viewpoint  has  distinct  advantages  which  when  applied,  are  comple- 
mentary. These  viewpoints  are  the  Statistical  Mechanics,  the  Structural, 

20 
and  the  Phenomeno log leal  approaches  to  material  behavior. 

The  Statistical  Mechanics  approach  considers  the  material  as  an 

association  of  discrete  particles  held  together  by  bonds  of  high  energy 

content.  At  this  level  behavior  of  the  discrete  particles  is  described 

by  their  relative  positions  in  space,  their  velocities,  and  the  Interaction 

forces  between  them.  The  Structural  approach  considers  the  material  to 

be  continuous  but  non-homogeneous,  being  formed  of  elements  of  different 

properties,  distributed  randomly  throughout  the  material,  and  having 

finite  dimensions.  The  Fhenomenologlcal  approach  considers  the  observed 

or  macroscopic  material  behavior.  It  deals  with  the  observation  of  the 

relation  between  forces  and  resulting  states  of  motion  of  finite  bodies 

assumed  to  be  homogeneous*  Miechanical  behavior  of  the  material  is 

described  in  terms  of  the  relations  between  stresses  and  strains,  and 

their  derivatives.  This  approach  usually  involves  extensive  testing  and 

bases  the  material  behavior  on  observed  responses  obtained  from  the  test 

and  its  environment. 
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A.   Objective 

In  this  Investigation,  an  attempt  has  been  made  to  combine  the  best 
findings  of  each  approach.  The  primary  concern  is  with  an  understanding 
of  the  behavior  of  a  specific  structural  lightweight  aggregate  concrete 
and  normal  limestone  aggregate  concrete  as  the  behavior  is  influenced  by 
creep  alone.  In  an  attempt  to  clarify  the  existing  theories  for  the 
mechanism  of  creep  as  they  pertain  to  concretes  of  the  type  tested,  this 
study  has  endeavored  to  link  molecular  hypothesis  and  rheological  models 
with  macroscopic  observations  of  the  concrete  structures  under  investi- 
gation. In  this  endeavor  it  is  hoped  that  a  clear  understanding  of 
creep  behavior  of  concrete  as  a  function  of  time  and  stress  will  result. 

The  parameters  in  the  model  developed  are  material  constants 
related  to  the  rheological  structure  of  the  concrete  in  question.  If 
it  is  possible  to  determine  the  nature  of  the  changes  of  these  material 
constants  for  the  concrete  structure  with  respect  to  stress  independently 
of  aging  it  is  reasonable  to  expect  that  it  will  be  possible  to  predict 
creep  for  the  concrete  under  any  condition  of  stress  and  time. 

In  all  previous  methods  developed  to  explain  creep  it  has  not  been 
possible  to  make  a  general  prediction  of  creep  behavior  for  a  specific 
concrete  relative  to  time  and  stress  based  on  the  results  from  a  single 
test.  The  advantages  for  being  able  to  do  so  are  obvious.  By  linking 
molecular  hypothesis  to  rheological  description  of  the  concrete  it 
appears  possible  to  separate  aging  effects  from  stress  variation  effects. 
When  the  two  effects  are  separated  variation  in  each  may  be  handled 
separately  thereby  producing  the  desired  model. 
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B.   Scope 

This  investigation  has  three  objectives. 

1.  To  determine  the  difference  in  the  nature  of  creep  between  lightweight 
Solite  concrete  and  normal  limestone  concrete. 

2.  To  determine  the  influence  of  shape  without  siae  interaction  on 
concrete  creep. 

3.  To  investigate  the  possibility  of  establishing  a  model  to  represent 
the  creep  of  concrete  under  any  stress  history  from  a  single  test 
of  the  concrete  which  is  to  be  represented  by  the  model. 

In  establishing  the  general  creep  model  with  respect  to  stress  and 
time  based  on  the  observed  creep  response  from  a  test  series  under 
sustained  stress,  conventional  viscoelastic  '   an(i  statistical  mechanics 
approaches  are  employed  to  develop  the  coefficients  of  the  creep  equation 
for  the  model.  Following  the  conventional  viscoelastic  approach  the 
material  behavior  is  separated  into  type  of  response  as  follows: 

1.  Instantaneous  elastic  response 

2.  steady  state  irrecoverable  creep 

3.  transient  delayed  elastic  response 

4.  transient  viscous  response  for  early  stage  creep 

5.  transient  irrecoverable  creep  representing  long  term  responses. 

The  model  is  established  by  converting  the  last  four  responses  listed 
into  three  parameters  in  the  model: 

1.  The  rate  of  transition  of  the  gel  from  fluid  to  solid. 

2.  The  total  expected  change  in  the  structure  of  the  gel  as  measured  by 
its  fluidity. 

3.  The  purely  elastic  transient  action  of  the  gel  and  aggregate  system. 
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The  first  parameter  constitutes  the  rate  of  change  of  the  concrete 
structure  under  the  action  of  stress  and  due  to  normal  aging.  The  second 
parameter  represents  the  difference  in  the  nature  of  the  structure  from 
before  load  is  applied  to  some  final  state.  These  two  parameters  are 
involved  in  the  stability  of  the  concrete  structure  with  respect  to 
aging  and  stress. 

The  third  parameter  is  purely  elastic  in  nature  and  represents  the 
delayed  elasticity  term  of  the  model.  It  is  shown  that  this  term  may  be 
employed  to  determine  the  delayed  elastic  recovery  for  the  unloaded 
specimen. 

If  the  model  is  to  be  acceptable  it  iiust  be  capable  of  predicting 
creep  behavior  of  lightweight  and  normal  concretes  alike;  the  model 
developed  does  so.  The  model  developed  also  attempts  to  predict  creep 
behavior  irrespective  of  size  and  shape  variations.  The  variations  in 
the  material  constants  of  the  model  due  to  shape  appear  to  adequately 
explain  the  creep  mechanism  as  it  pertains  to  shape.  Size  was  not  a 
parameter  of  the  experiment.   However  the  anticipated  nature  of  the 
interaction  between  size  and  shape  is  developed  and  is  presented  in  the 
interpretation  of  the  results.   Since  mix  variations  and  environment  were 
not  parameters  in  this  experiment  and  in  view  of  a  lack  of  sufficient 
correlation  between  size  and  creep,  the  scope  of  this  experiment  precludes 
the  formulation  of  a  more  generalized  model. 


CHAPTER  II 

REVIEW  OF  THE  CREEP  BEHAVIOR  OF  CONCRETE 

Many  variables  influence  the  magnitude  of  creep.  This  test  program 
was  limited  to  only  a  few  among  which  are: 
Coarse  aggregate  Ingredients 
Mix  proportions 
Shape  of  the  specimen 

A,  Coarse  Aggregate  Ingredients 

In  general  it  is  believed  that  hard  dense  aggregates,  with  low 
absorption  and  a  high  elastic  modulus  produce  a  concrete  with  low  creep 
tendencies.  Troxell  and  Davis  Indicate  that  particle  shape,  surface 
texture,  pore  structure  and  unit  weight  may  also  Influence  concrete  creep, 
These  influences  are  all  involved  in  the  differences  between  the  two 
aggregates  used  in  this  study.  An  individual  evaluation  of  each 
influencing  factor  is  difficult  if  not  impossible.  However,  it  may  be 
possible  to  speculate  on  the  more  prominent  characteristics  which  are 
considered  to  be  the  modulus  of  elasticity  and  the  permeability  of  the 

aggregate.  Particle  shape  and  surface  texture  are  not  considered  to 

15 
be  major  factors  in  creep  behavior  according  to  Best. 

B.  Mix  Proportions 

Cement  paste  content  and  water-cement  ratio  appear  to  be  intimately 
involved  in  creep  activity  of  concrete.  Many  authors  have  testified  to 
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the  direct  relationship  between  paste  content  and  creep.    This  study 

was  not  intended  to  resolve  the  questions  concerning  the  nature  of  the 

involvement  of  cement,  water  or  paste.  Therefore,  in  an  attempt  to 

minimize  complications  in  the  behavior  of  concrete  due  to  variations  in 

these  ingredients  the  paste  contents  were  made  equal  for  the  concrete 

mixes . 

The  water-cement  ratios  for  the  two  mixes  could  not  be  so  easily 

equated.   No  problem  resulted  from  the  selection  of  a  given  water-cement 

ratio  for  the  limestone  concrete.  However,  a  water-cement  ratio  for  the 

lightweight  aggregate  concrete  may  not  be  nstablished  since  a  portion  of 

the  free  water  contained  in  the  aggregate  used  must  be  considered  to  add 

to  the  mix  water  for  the  concrete.   By  preliminary  experimentation  the 

strength  of  the  lightweight  concrete  was  controlled  by  selection  of  an 

effective  mix  water  content  to  establish  near  equal  strengths  for  the 

two  structural  concretes. 

C.   Shape  of  Specimen 

Very  little  work  has  been  done  to  determine  the  effect  of  shape  on 
creep.  Several  in  vestigations  have  been  reported  involving  the  influence 
of  size  of  specimen  on  creep.  These  reports  indicate  a  decrease  in 
creep  with  increase  in  size  of  specimen.^'   '     The  effect  is  usually 
explained  as  a  result  of  the  increased  seepage  path  to  the  free  surface. 
This  apparent  interaction  of  stress  dependent  seepage  and  shrinkage  which 
is  related  to  creep  bears  directly  on  the  shape  effect  of  the  concrete 
specimen.  Thus  by  varying  the  free  surface  area  and  maintaining  a 
constant  cross-sectional  area  it  is  possible  to  evaluate  the  stress 
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dependent  seepage  effect  on  creep  as  distinguished  from  pure  shrinkage 
for  these  two  concrete  types. 

D.  Nature  of  Creep  in  Concrete 

The  deformation  characteristics  of  concrete  are  related  to  its 

11  13 
internal  structure.   '    The  mortar  paste  constitutes  the  active  part 

of  the  concrete  in  binding  coarse  aggregate  to  form  a  stone-like  material. 

It  is  composed  of  cement,  water,  interspersed  graded  fine  aggregate,  and 

air.  The  relative  proportions  of  these  Ingredients  vary  from  one  concrete 

mix  to  another.  This  variation  is  the  essential  factor  which  results  in 

all  behavioral  variations  of  concrete. 

The  aggregate  particles  are  bound  together  by  the  cement  paste 

matrix.  The  aggregate  plays  a  less  active  role  in  the  creep  response 

of  the  concrete  than  the  cement  paste. 

At  any  stage  of  hydration,  the  structure  of  the  cement  paste 
consists  of  unhydrated  cement  grains,  hydrated  cement,  water,  and  air. 
The  cement  paste  structure  is  called  cement  gel  which  has  a  more  rigid 
structure  with  aging  and  a  less  rigid  structure  under  increasing  load. 
It  Is  the  action  of  the  gel  which  gives  concrete  its  viscoelastic 
nature.  Initially  as  the  cellular  structure  of  the  hydrated  cement 
is  developing,  the  colloidal-liquid  phase  predominates.  If  a  load  is 
applied  at  this  early  stage  the  cement  gel  behaves  as  a  liquid  in  that 
it  flows  under  any  nonisotropic  stress,  however  small.   The  degree  to 
which  the  liquid-colloidal  phase  exists  is  manifested  in  the  amount  of 
creep  which  results.  As  the  gel  matures,  due  to  continued  hydration, 
the  solid  cellular  structure  becomes  dominant  over  the  liquid  phase. 


An  applied  stress  at  this  time  results  in  less  creep  since  the  gel  would 
behave  more  like  a  solid,  exhibiting  elastic  tendencies. 

The  nature  of  the  action  of  the  aggregate  dispersed  throughout  the 

2 
cement  gel  has  been  termed  passive.   If,  as  stated  above,  it  is  the 

cement  gel  which  creeps,  then  by  including  solid  aggregate  particles 

into  the  gel  its  resistance  to  creep  should  be  increased  according  to 

3 
Reiner.   This  condition  further  complicates  the  problem  of  determining 

the  action  of  the  coarse  aggregate.  As  it  was  explained  in  connection 

with  the  coarse  aggregate  effects  discussed  in  section  A  of  this  chapter, 

this  additional  volume  concentration  effect  cannot  be  evaluated  for  the 

lightweight  aggregate  because  it  is  not  known  how  the  moisture  contained 

in  the  lightweight  aggregate  influences  its  volume  concentration  in  the 

concrete . 

Therefore,  like  most  creep  findings  reported  in  the  literature, 

the  findings  of  this  report  may  only  be  applied  to  the  specific 

concretes  being  investigated.  The  model  constants  are  applicable  only 

to  the  range  of  specimen  sizes  used  and  under  similar  ambient  conditions 

of  humidity  and  temperature. 


CHAPTER  III 

THE  VISCOELASTIC  APPROACH 

A.  Basic  Mathematical  Models 

The  most  successful  simple  phenomenologlcal  £orm  of  the  data  from 
a  creep  study  follows  Ross's  parametric  equation  relating  strain  and 
time.  Most  strain  data  may  first  be  prepared  in  this  form  from  which 
creep,  e  ,  may  be  determined  at  any  time  t  from 

at 


b  +  t 


(1) 


where  a  is  a  dimensionless  constant  and  b  is  a  constant  with  the  units 
of  time. 

The  most  general  phenomenologlcal  form  for  creep  data  is  obtained 
from  rheologlcal  theory.  Rheology  is  the  science  of  flow  of  materials; 
it  offers  a  transition  between  classical  theory  of  elasticity  and  classical 
hydrodynamics  and  has  been  proposed  as  a  theory  of  the  general  behavior 
of  materials  on  the  assumption  that  every  real  material  must  be  supposed 
to  possess  all  basic  deformatlonal  properties  in  varying  proportions. 

Strain  relations  may  be  developed  as  viscoelastic  strains  employing  the 

2  17 
conventional  theory  of  rheology.  *     Using  the  data  from  test  series  A 

a  creep  model  with  four  linear  elements  was  developed  which  was  similar 

to  Freudenthal ' s  model  in  the  simplified  form.  The  detailed  model  used 

is  shown  in  Figure  1.  The  model  consists  of  the  following  elements  in 

series . 
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LINEARIZED     RHE0L06ICAL   MODEL  FOR   CONCRETE 

FIGURE  1 
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1.  A  Maxwell  element  with  linear  spring  and  dashpot  in  series  representing 
the  elastic  strain  and  the  long  time  irrecoverable  creep  resulting  from 
inelastic  strain. 

2.  A  Kelvin  element  with  linear  spring  and  dashpot  in  parallel  representing 
the  short  time  viscoelastic  effects.   This  Influence  was  attributed  to 
consolidation  due  to  seepage  of  pore  water. 

3.  A  Kelvin  element  with  linear  spring  and  dashpot  in  parallel  representing 
intermediate  time  viscoelastic  effects.   This  influence  was  attributed, 
after  Freudenthal,  to  retarded  elasticity  or  recoverable  creep. 

4.  A  Kelvin  element  with  linear  spring  and  dashpot  in  parallel  representing 
long  time  viscoelastic  effects.   This  Influence  was  attributed  to 
destruction  of  the  gel  structure  under  stress. 

The  equation^'  representing  the  simple  rheologlcal  model  for  constant  f  is 

e  =  |^lt_  +,     I-  (l-e-^nVj  ^^^ 

^       K  n  =  1   n 

The  first  two  terms  represent  the  Maxwell  element.  The  three  terms  in 
the  summation  represent  the  Kelvin  elements.   All  terms  in  equation  (2) 
are  defined  In  the  list  of  symbols. 

Most  creep  curves  may  be  satisfied  by  simply  extending  the  number 
of  terms  in  the  summation  to  provide  a  transition  curve  which  fits  the 
data  between  the  elastic  response  (first  term)  and  the  final  or  steady 
state  response  (second  term). 

B.   The  Nature  of  the  Rheo lop leal  Model 

Rheologlcal  equations  define  ideal  bodies  which  serve  as  models  of 
comparison  in  the  analysis  of  material  behavior.  The  rheologlcal 
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varlables  are  stress,  deformation  and  time;  the  rheologlcal  parameters 
are  viscosity  and  elasticity.  The  parameters  in  the  rheologlcal  equations 
characterize  the  material  behavior;  they  may  be  constants  £rom  which  we 
obtain  linear  vlscoelastic  equations  or  they  may  be  functions  of  time  or 
stress  from  which  we  obtain  non-linear  vlscoelastic  equations. 

Many  models  have  been  proposed  for  the  creep  behavior  of  concrete. 
However,  the  many  factors  influencing  the  gel  structure  seem  to  preclude 
the  development  of  a  single  model  for  all  conditions.  The  nature  of 
hydration  in  the  gel  is  believed  to  contribute  the  major  source  of 
variability  in  the  creep  behavior.  It  is  the  gel  which  supplies  the 
vlscoelastic  nature  of  concrete  because  it  is  a  varying  form  of  matter 
lying  between  the  solid  and  fluid  physical  states.  The  gel  in  concrete, 
which  is  formed  Initially  from  a  fluid,  is  distinguished  by  a  change  in 
Its  mechanical  properties,  by  a  transition  into  a  solid  state  having  high 

viscosity,  elasticity,  and  limiting  values  of  stress  related  to  strain. 

5 
Some  authors  have  experienced  completely  solid  state  responses  (no 

viscous  creep)  from  concretes  which  have  been  fully  dried  and  cured. 

Concrete  undergoes  the  transition  from  fluid  to  solid  over  a  period 

of  several  years.  After  a  few  hours,  however,  the  solid  phase  cellular 

structure  may  have  developed  sufficient  elastic  strength  to  withstand  a 

test  load.  At  that  time  the  concrete  is  put  into  service  and  loads  are 

imparted  to  it.  It  is  clear  that  the  degree  of  fluidity  of  the  gel  at 

the  time  of  loading  will  greatly  influence  the  elastic  and  creep  properties 

of  the  concrete.  Rheologlcal  models  are  unable  to  cope  with  the  many 

variables  which  influence  these  phase  changes  of  the  gel  not  to  mention 

the  added  variation  of  aggregate  interaction  with  the  gel.  It  would 
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appear  that  a  single  model,  no  matter  how  general  it  may  be,  is  not  the 
basic  answer  to  prediction  of  creep  in  concrete;  rather,  a  general 
procedure  for  predicting  creep  behavior  is  more  in  order.   Such  a  procedure 
will  provide  for  the  development  of  an  equation  which  incorporates  the 
significant  parametiers  for  the  concrete  being  investigated  into  a  specific 
rheo logical  model. 

All  the  influences  on  creep  may  be  categorized  by  their  influence 
on  three  factors: 

1.  The  rate  of  transition  of  the  gel  from  fluid  to  solid. 

2.  The  purely  elastic  action  of  the  gel  and  aggregate  system. 

3.  The  total  change  in  the  structure  of  the  gel. 

If  these  three  factors  are  used  as  parameters  for  the  rheo logical 
model  then  it  would  seem  that  the  creep  behavior  of  every  concrete  may 
be  interpreted . 

For  a  more  detailed  study  of  the  rheological  approach  refer  to 
items  2,  4,  and  17  in  the  bibliography.   A  supplementary  study  to  this 
report,  of  rheology,  and  the  basis  for  development  and  behavior  of 
the  model  presented  in  figure  1  is  found  in  reference  22. 


CHAPTER  IV 
THE  STATISTICAL  MECHANICS  APPROACH 

A.  The  Rate  Process  Theory 

6 
Glasstone,  Laidler  and  Eyring  extended  classical  quantum  mechanics 

theory  to  include  statistical  treatment  of  reaction  rates.  Since  all 

matter  is  composed  of  molecules  or  assemblies  of  molecules  which  vibrate 

about  some  equilibrium  position,  when  the  quantum  state  (or  energy  level) 

is  changed  new  equilibrium  positions  are  obtained.  Eyring,  et.  al .  employed 

a  concept  of  potential  energy  surfaces  to  describe  the  conversion  from 

vibrational  energy  to  translational  kinetic  energy  and  vice-versa.  The 

essential  requirements  in  the  development  are  that: 

1.  Energy  is  conserved,  following  the  first  law  of  thermodynamics. 

2.  The  number  of  assemblies  is  constant. 

3.  There  exists  definite  energy  levels  (second  law  of  thermodynamics). 

4.  All  possible  energy  levels  (or  quantum  states)  for  the  entire  system 
have  equal  probability. 

Since  there  is  a  given  probability  that  any  molecule  or  assembly 

shall  have  a  free  energy  in  any  quantum  state  which  is  a  function  of  the 

entropy  of  the  assembly,  the  total  probability  of  the  occurrence  of  an 

assembly  with  given  entropy  is  proportional  to  the  sum  of  all  the  energy 

terras  for  the  assembly.  The  sum  is  called  the  partition  function  of  the 

assembly  for  a  given  volume  of  matter.  The  partition  function  may  include 

energy  terms  for  nuclear,  electrical,  vibrational,  rotational,  and 

translational  energy. 
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In  evaluating  the  reaction  rate  of  an  assembly  we  are  concerned  about 

the  velocity  at  which  an  activated  assembly  travels  over  the  potential 

9 
energy  barrier  (Figure  2)  thus  passing  from  one  equilibrium  state  to 

another.  The  magnitude  of  the  energy  barrier  is  equal  to  the  work  or 

thermal  energy  which  must  be  applied  to  the  assembly  in  order  to  activate 

it.  The  assembly  is  considered  to  be  activated  when  it  is  at  the  level 

of  energy  equal  to  the  energy  barrier  level.  The  net  rate  at  which  the 

reaction  occurs  is  determined  by  the  average  velocity  of  the  activated 

assemblies  passing  over  the  top  of  the  barrier. 

Eyring  et  al.  have  evaluated  the  partition  functions  involved  in 

the  thermodynamics  of  reaction  rates  to  develop  an  equation  for  the 

specific  reaction  rate.  The  specific  reaction  rate  defines  the 

frequency,  v,  with  which  an  activated  assembly  crosses  the  barrier  and 

is  displaced  a  distance  A  (Figure  2)  from  one  equilibrium  position  to 

another. 


e  ■  ^°'"  (3) 


kT  .   "^o 
h 


in  which: 

T  "  absolute  temperature 
k  a  Boltzmann's  constant 
h  =»  Planck's  constant 

E  =  the  free  energy  of  activation;  for  a  rheological  assembly  it  is 
defined  by  E^  =  V  -  T  •  S  where  V  is  the  total  energy  of  the 
assembly  and  S  is  the  entropy  (a  measure  of  the  disorder  of  the 
assembly)  accompanying  the  activation  process. 
In  the  presence  of  an  applied  stress  the  energy  barrier  is  modified 
(Figure  2).  The  effect  of  the  stress  is  to  reduce  the  height  of  the 
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energy  barrier  in  the  direction  of  the  stress  thereby  permitting  a  flow 
assembly  with  given  free  energy  level  to  pass  freely  over  the  barrier. 
Since  there  exists  an  equal  probability  for  the  existence  of  all  free 
energy  levels,  a  greater  stress  will  result  in  activation  of  a  greater 
portion  of  the  population  of  flow  assemblies. 

If  the  potential  energy  surface  is  considered  symmetrical,  then  the 
energy  level  required  for  flow  opposite  to  the  applied  stress  is 
increased  an  amount  equal  to  the  decrease  in  the  energy  level  in  the 
direction  of  the  applied  stress.  That  is,  activation  occurs  midway 
between  equilibrium  positions. 

The  passage  of  an  activated  assembly  over  the  potential  energy 
barrier  represents  the  jump  of  the  assembly  from  one  equilibrium 

position  to  the  next.  Let  ^  be  the  distance  between  two  equilibrium 

2 
positions.  The  applied  force  acting  on  an  assembly  is  f/D  where  D 

is  a  structure  factor  which  is  equal  to  the  number  of  flow  assemblies 

per  unit  length  and  f  is  the  applied  stress.  Hence,  the  energy  that  the 

assembly  acquires  in  advancing  to  the  activated  state  is 

o  2 

Af/D  X  1/2    ,  that  is,  1/2  f A  /D  . 

A  creep  condition  may  be  analyzed  by  evaluating  the  coefficients  for 

the  specific  reaction  rate  of  the  concrete  specimen  loaded  under  sustained 

stress.  The  frequency  of  activation  for  movement  of  the  flow  assemblies 

in  the  direction  of  the  stress,  that  is,  for  flow  in  the  forward  direction, 

is 

Vf  =  -IT  exp     '_° 2Df  <3a) 
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and  the  frequency  for  movement  opposed  to  the  stress,  that  is,  for  flow 
in  the  backwards  direction,  is 


kT 


■-(Ko.fs^) 


(3b) 


kT 

The  net  frequency  of  movement  of  the  flow  assemblies  combining  motions 
in  the  forward  and  backward  directions  becomes 

kT     l""^  I       f  A 
Vf  -  ^b  =  2  "h  exp  [  kT  Jsinh  (2^x0^^   '  ^^""^ 

The  creep  rate  may  now  be  related  to  the  frequency  of  movement  of  flow 
assemblies.  Let 

L  =  the  axial  component  of  displacement  due  to  the  movement  of 
an  assembly  to  a  new  equilibrium  position,  and  recall  that 
D  =»  a  structure  factor  which  is  equal  to  the  number  of  flow 
assemblies  per  unit  length. 
The  creep  rate  is 

e=  DL  (vf  -  Vb) 


kT  [""^o"!  fA 

€   =  2  DL~    exp    |_kTjsinh  (2^^)  (^) 

which  may  be  written 

e  =  A  sinh  B  f       .  (4a) 

The  coefficients,  A  and  B  are  evaluated  based  on  the  findings  of 
the  rheological  investigation  of  the  parametric  equations  obtained  from 
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the  observed  responses.  Using  the  simple  viscoelastic  model  of  Figure  1 
described  in  Chapter  III  with  the  rate  process  solution,  the  data  may  be 
further  evaluated  for  a  better  understanding  of  the  creep  behavior. 

The  rate  process  equation  (4a)  relates  creep  rate  to  stress.  The 
viscoelastic  equation  (2)  for  the  simple  rheological  model  of  Figure  1 
related  creep  to  time  for  a  given  sustained  stress.  The  two  applications 
are  therefore  complemttntary  and  may  facilitate  separation  of  influences 
on  creep  due  to  aging  and  stress  variations. 


CHAPTER  V 

EXPERIMENTAL  PROCEDURE 

A.  Description  .of  the  Test 

It  is  generally  considered  that  creep  and  shrinkage  are  Intimately 
related.  Their  true  Interaction,  however,  is  not  completely  understood. 
Shrinkage  has  been  defined  as  the  change  in  length  of  concrete  members 
without  the  influence  of  applied  stresses.  Creep  is  defined  as  the 
change  in  length  of  concrete  members  under  the  influence  of  applied 
stresses.  These  definitions  provide  only  a  superficial  distinction 
between  creep  and  shrinkage.  In  fact,  differential  shrinkage  in 

concrete  is  known  to  result  in  stress  gradients  causing  compressive 

21 
stresses  in  the  interior  of  the  concrete.    Therefore  any  contribution 

to  the  shortening  of  the  member  by  these  stresses  should  be  classified 

as  creep  according  to  the  definition.  On  the  other  hand,  moisture  forced 

to  the  vicinity  of  the  surface  of  a  creep  specimen  by  applied  stresses 

may  be  removed  from  the  specimen  by  the  same  phenomena  that  cause  shrinkage. 

In  spite  of  this  awareness,  it  is  traditional  to  separate  the  two  effects, 

shrinkage  and  creep,  according  to  the  definitions  presented  above.  This 

separation  is  desirable  in  this  investigation  since  the  basic  theories 

being  employed  are  contingent  upon  a  knowledge  of  the  applied  forces. 

Consequently  creep  strains  in  this  investigation  were  obtained  by  subtraction 

of  shrinkage  strains  obtained  from  unstressed  specimens,  from  the  total 

strains  obtained  from  creep  specimens. 

-20- 
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The  purpose  of  this  Investigation,  concerning  the  comparative  creep 
behavior  of  a  typical  Florida  limestone  aggregate  concrete  and  a  Florida 
lightweight  aggregate  concrete  (Florida  Solite  -  an  expanded  clay  product) 
was  to  determine  some  structurally  Important  aspects  of  the  mechanism  of 
creep  in  concrete.  Three  test  series  investigated  creep.  Series  A 
investigated  creep  under  sustained  stress.  Test  series  B  and  C  Investigated 
creep  under  different  decreasing  stress  situations.  A  fourth  test  series 
investigated  shrinkage. 

Shrinkage  of  the  test  specimens  is  illustrated  in  Figures  12,  13,  and 
14.  Total  strains  for  creep  series  A  are  Illustrated  in  Figures  IS,  16, 
and  17.  Both  sets  of  strain  data  display  seasonal  influences  beyond  the 
initial  stage  of  their  responses.  Shrinkage  and  creep  specimens  were 
maintained  in  the  laboratory  environment  where  periodic  ambient  readings 
were  recorded  of  humidity  (Figure  6)  and  temperature  (Figure  7). 

B.  Materials 

In  an  attempt  to  minimize  effects  from  mix  variations,  strength, 
variations  and  environment  variations  between  pours,  a  single  pour  was 
made  for  each  concrete  type.  Table  I  contains  proportions  for  the  two 
mixes  used.  Equal  paste  contents  as  nearly  as  could  be  determined  were 
selected  to  minimize  the  effects  of  a  difference  between  the  concrete 
gels  influencing  creep.  As  a  result  of  the  desire  to  obtain  near  equal 
strength  for  the  two  concrete  types,  the  water  contents  were  selected  to 
compensate  for  the  aggregate  influence  on  strength.  Figure  3  Illustrates 
the  concrete  strength  gain  for  both  concretes.  Each  point  represents  the 
average  of  four  6  x  12  inch  test  cylinders. 
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— -The  concrete  used  in  the  specimen  contained  type  III,  high  early 
strength  cement  as  would  be  used  in  most  prestressed  concrete  bridge 
construction.  The  coarse  aggregates  had  maximum  sizes  o£  3/4".  The  fine 
aggregate  used  in  all  specimens  was  Interlachen  sand  normally  used  for 
construction  in  the  north  Florida  area. 

C.  Description  of  Specimens 

Special  shapes  were  selected  to  test  an  hypothesis  about  the  shape 
factor  effect  on  the  rate  of  creep.  Three  shapes  with  identical  cross- 
sectional  areas  equal  to  28.27  square  inches  but  with  unlike  perimeters 
were  cast.  The  specimens  obtained  contained  equal  volumes  but  varying 
volume-to-surface  area  ratios  (Table  II).   The  shapes  selected,  which 
had  cross,  rectangular,  and  circular  cross-sections,  are  shown  in  Figure  4. 
The  selection  of  equal  volumes  was  made  in  order  to  evaluate  shape  effect 
by  blocking  out  size  effect. 

The  flow  path  from  any  point  at  the  interior  of  the  specimen  to  the 
surface  is  measured  by  the  average  weighted  distance  over  which  a  particle 
of  moisture  must  travel  before  reaching  the  surface  of  the  specimen.  The 
average  weighted  distances  for  the  specimen  may  be  calculated  by  summing 
the  products  of  the  incremental  areas  of  the  specimen  cross-section  and 
the  least  distance  from  these  incremental  areas  to  the  surface  of  the 
specimen,  and  dividing  the  sum  by  the  area  of  the  specimen.  The  average 
weighted  distances  are  given  in  Table  II  for  each  shape  of  specimen. 
Although  perimeter-to-cross-sectlon  area  ratio  is  frequently  used  as  a 
measure  of  shape  effect   it  is  convenient  and  perhaps  more  accurate  to 
use  the  weighted  distance  as  the  significant  measure  of  shape  effect. 
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When  this  substitution  is  made,  it  becomes  immediately  evident  that  shape 
and  size  effects  interact  since  length  of  flow  path  is  the  primary 
variable  governing  their  influences  on  creep.  In  this  regard  the  specimen 
behavior  in  this  investigation  may  be  considered  to  give  some  indication 
of  size  effect  even  though  the  investigation  seemingly  blocks  out  size 
effects. 

All  specimens  were  thirty  inches  long.  This  provided  space  for 
two  gage  lengths  of  ten  inches  over  which  measurements  were  made  plus 
six  inches  on  each  end.  The  six  specimens  comprising  a  creep  test  series 
were  strung  together  end  to  end  with  a  transition  plate  between  each  specimen. 
The  ends  of  all  specimens  were  capped  with  standard  capping  compound 
consisting  of  flyash  and  sulfur.  Each  specimen  was  cast  with  a  metal 
conduit  along  its  axis  through  which  the  stress  rod  was  placed  for 
loading  of  the  specimen.  A  ten  inch  concrete  bearing  clock  and  steel 
plate  were  used  at  the  ends  of  the  strung  specimens  to  accept  and 
distribute  the  applied  test  load. 

Shrinkage  specimens  were  made  identical  to  creep  specimens.  Their 
end  faces  were  sealed  with  a  layer  of  wax  and  metal  foil  which  resulted 
in  having  their  exposed  surface  identical  to  the  exposed  surface  of  the 
creep  specimens. 

Shrinkage  strains  were  taken  on  individual  unstressed  specimens. 
Creep  stains  were  taken  on  specimens  under  identical  stress  conditions 
in  each  test  series.  One  specimen  of  each  shape  for  each  of  the  two 
concrete  types  was  used  in  each  creep  teat  series  as  shown  in  Figure  5. 
Table  II  summarizes  the  test  conditions  for  each  specimen.  Specimens  In 
test  series  A  were  maintained  under  constant  sustained  stress  of  914  -  IQ  psi. 
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Stresses  in  specimen  for  test  series  B  and  C  were  permitted  to  decrease 
as  shown  in  Figure  8.  The  decrease  in  stress  associated  with  series  C 
resulted  from  normal  prestress  loss  due  to  creep  and  shrinkage.  Test 
series  B  contained  filler  blocks  of  concrete  and  high  strength  plaster 
between  the  specimens  which  permitted  a  greater  initial  stress  decay. 

D.  Instrumentation 

All  specimen  were  instrumented  with  standard  Whittemore  gage  points 
spaced  at  ten  inches  as  indicated  in  Figures  4  and  5.  Shrinkage  and 
creep  specimens  were  maintained  in  laboratory  environment  where  periodic 
concrete  strain  readings  were  recorded  along  with  ambient  readings  of 
humidity  (Figure  6)  and  temperature  (Figure  7).  Strain  adjustments 
were  made  for  temperature  variation  by  assuming  a  thermal  coefficient  of 
0.0000055  inches  per  inch  per  degree  Fahrenheit.  However,  no  adjuetmente 
were  made  for  humidity  variations  because  there  is  no  way  of  determining 
the  Integration  effects  of  humidity  on  concrete  behavior. 

Load  readings  on  the  creep  specimen  were  obtained  from  SR-4  strain 
gages  on  the  stress  rods.  Loads  were  applied  with  hydraulic  Jacks. 
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CHAPTER  VI 


DEVELOPMENT  OF  A  MODEL  FOR  THE  CREEP  MECHANISM 


A.  Separation  of  Creep  Components 

There  are  two  fundamental  types  of  viscoelastic  creep:, 

1.  Delayed  elasticity  --  a  completely  recoverable  strain  phenomenon. 

2.  Viscous  creep  ~  an  irrecoverable  strain  which  may  contain  transient 
and  steady  state  components. 

Creep  has  been  associated  with  rearrangement  of  molecules  due  to 
thermal  movement.  According  to  quantum  theory,  the  movement  must  result 
in  a  different  and  stable  equilibrium  configuration  of  molecules.  Since 
fluids  contain  many  vacant  sites  in  their  molecular  structure,  large 
voids  (holes)  exist  into  which  the  migrating  molecules  may  relocate.  Any 
motion  of  one  molecular  assembly  may,  of  course,  result  in  relocation  of 
neighboring  assemblies.  The  rearrangement  of  these  assemblies  is  always 
in  the  general  direction  of  the  applied  stress  but,  the  motion  may  have 
components  in  any  direction  thereby  resulting  in  shear  distortion.  These 
shear  distortions  are  viewed  as  macroscopic  creep.  When  the  applied  stress 
is  removed  some  of  these  migrated  assemblies  return  to  their  initial 
equilibrium  states.  The  resulting  strain  recovery  is  viewed  as  delayed 
elasticity.  It  is  necessary  to  separate  the  two  fundamental  types  of 
creep,  the  irrecoverable  from  the  recoverable  creep,  in  the  development 
of  the  model  for  a  creep  mechanism  because  under  decreasing  stresses  the 
elastically  recoverable  component  will  be  active. 
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(a)  Elastic  Creep 

Q 

Elastic  creep  has  been  described  by  Orowan  as  resulting  from  an 
applied  stress  which  alters  the  statistical  equilibrium  position  £or  each 
assembly.  A  time  dependent  strain  results  which  continues  as  the  assemblies 
f  seek  new  equilibrium  positions  within  their  stressed  configuration.   If 

the  force  is  suddenly  removed,  the  instantaneous  elastic  response  returns 
the  molecular  structure  to  its  former  equilibrium  configuration.  The 
assemblies  now  seek  their  original  equilibrium  positions  and  recoverable 
strains  are  observed.  It  is  tacitly  assumed  that  no  deterioration  or 
damage  occurs  to  the  molecular  structure  during  the  load  cycle  and  no 
additional  thermal  energy  sources  are  introduced  into  the  system.  The 
strain  recovery  will  therefore  be  complete. .by  definition  of  elasticity. 

(b)  Viscous  Creep 

( 

During  the  transition  period,  while  fluidity  characterizes  the  nature 

of  the  gel  matrix,  the  flow  assemblies  are  more  mobile.  The  crystalline 

structure  is  not  sufficiently  rigid  to  behave  completely  elastically  and 

therefore,  once  it  has  moved  to  a  new  equilibrium  position  the  flow  unit 

is  completely  stable  in  this  position.  Upon  removal  of  the  load  there  is 

no  tendency  for  the  flow  unit  to  return  to  Its  original  position, 

B.  The  Generalized  Model  for  the  Creep  Mechanism 
I  The  development  of  a  model  which  Is  to  represent  the  creep  behavior 

of  concrete  must  contain  the  two  elements,  delayed  elasticity  and  viscous 

14 
creep.  The  basic  strain  equation  must  then  be  of  the  form 

€=  E  +  q  (f.  t)vis.  +  *1  <^'  '^^el.    .  ^^^ 

where  the  terms  represent  the  elastic  strain,  the  viscous  component  of 
creep,  and  the  elastic  conqponent  of  creep  respectively. 
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Since  the  last  term  has  an  elastic  nature,  a  Kelvin  element  should  be 
sufficient  to  represent  this  component.  In  its  simplest  form  the  term  may 
have  constant  coefficients  to  represent  average  effects  of  delayed  elasticity 
in  the  concrete.  It  represents  the  basic  elastic  structure  of  the  gel. 

The  second  term  contains  all  phenomena  which  result  in  permanent  changes 
in  the  rheological  structure.   Such  phenomena  are  for  example: 

1.  Hydration  effects. 

2.  Irrecoverable  strain  effects  resulting  from  structural  deterioration 
due  to  stress. 

3.  Viscous  flow. 

Each  of  the  two  creep  terms  in  the  model  must  contain  coefficients 
which  relate  the  structural  changes  of  the  concrete  due  to  the  various 
influences  to  the  two  major  variables,  aging  and  stress.  The  objective  is 
to  establish  coefficients  which  include  as  many  influencing  factors  as 
possible  so  as  to  make  the  model  applicable  in  a  greater  number  of  situations. 
For  example,  if  it  would  be  possible  to  incorporate  such  factors  as  size 
and  shape,  environment  conditions,  mix  proportions,  and  physical  characteristics 
of  the  mix  ingredients  into -the  model,  then  the  enigma  of  creep  behavior  would 
truly  be  considered  resolved.  It  should  be  emphasized  at  this  point  that 
•  this  study  will  not  attempt  to  develop  such  a  general  form  of  the  model 
for  the  creep  mechanism.  It  is  hoped,  however,  that  the  ground  work  will 
be  adequately  accomplished  toward  such  a  development.  Obviously  to  achieve 
such  a  result  will  involve  a  study  of  much  greater  depth  and  scope  than 
this  investigation  encompasses. 

Equation  (5)  may  be  revised  to  include  the  delayed  elastic  component 

referred  to  and  is  therefore  written  in  the  form: 

c  f  -0eEe'=. 

C=  |-  +  0  (f,  t)  ft  + 1-  (1  +  e      )    .  (5a) 


-36- 


(a)  Coefficients  of  Structural  Stability 

3 
Reiner  has  introduced  a  coefficient  (X)  called  "The  Coefficient 

of  Structural  Stability"  i^lch  relates  viscosity,  or  more  accurately  Its 

reciprocal,  fluidity, to  stress.  The  scope  of  this  coefficient  may  be 

broadened  to  include  structural  variations  in  fluidity  due  to  aging  also. 

As  described  in  Appendix  B  let 

where    0     represents  the  fluidity  at  any  time, 
000    represents  the  ultimate  fluidity, 

f     is  the  normal  stress. 

2 

The  partial  derivatives  are  taken  with  respect  to  stress  squared  (f  ) 

and  time  (t).  Reiner  has  shown  that  fluidity  is  Influenced  by  the  square 

2 
of  the  shear  stress  (f  /4).  The  negative  sign  is  required  because  an  increase 

in  time  results  in  a  decrease  in  fluidity. 

In  order  to  understand  the  Structural  Stability  Coefficient  consider 

the  extreme  cases  of  X  »  0  and  X  »  °°  ,  When  X  «  0,  0  "  0„      which  means 

that  the  viscosity  of  the  material  Is  at  its  steady  state  value.  The 

material  is  fully,  aged  and, under  applied  stress,  it  strains  at  constant  rate. 

When  X  <°>  00  ,  the  partial  derivatives  of  fluidity  with  respect  to  time  and 

stress  are  equal  to  zero.  This  Implies  that  the  viscosity  of  the  material 

is  not  changed  due  to  applied  stress  or  due  to  aging.  This  latter  case  Is 

obviously  not  possible  with  concrete  since  aging,  through  hydration,  will 

always  produce  some  finite  value  for  the  change  in  viscosity  with  time. 

Therefore  the  smaller  values  for  the  structural  stability  coefficient,  X, 

indicate  greater  structural  stability  of  the  material.  The  material  Is  less 

subject  to  creep,  since  0  is  smaller  because  of  aging  or  stress  Infludncea. 
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Either  factor,  stress  or  aging,  may  have  a  greater  influence  or  a 
lesser  influence  on  structural  changes.  It  is  therefore  necessary  to 
separate  the  aging  and  stress  components  of  the  coefficient  of  structural 
stability.  Accordingly  components  for  aging,  c  ,  and  stress  change,  C£, 
have  been  developed  in  Appendix  B. 

It  can  be  shown  (Appendix  B)  that  the  coefficient  of  structural 


stability  will  influence  the  fluidity  in  the  following  way: 

2   ,2. 


,)  exp  [-  (ct  +  Of  (fo  -  f))]  , 


(7) 


where 

000  =  steady  state  or  ultimate  fluidity 

0o  =  initial  fluidity 

fo  =  initial  stress 

f  =  stress  at  any  time  (t) 

The  coefficient  of  structural  stability  contains  two  components: 

1,  Cf  represents  the  extent  of  influence  of  stress  change  on  the  structural 
stability. 

2.  c  is  Independent  of  stress  change  and  represents  the  natural  influences, 
such  as  hydration,  on  the  structural  stability. 

Substituting  equation  (7)  into  equation  (5a),  the  strain  may  now  be 
written 

c  4  +[«^co  +  (K-K)  exp  (-  (c,  +  c^  (C-f')))]  ft  +  i;[i-e-'*^^«*]/8) 

This  equation  represents  a  general  function  for  creep  following  the  mechanical 
model  of  Figure  9.  It  becomes  necessary  to  evaluate  the  coefficients  an 
parameters  of  the  second  and  third  elements  in  Figure  9.  The  parameters  appear 
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in  the  form  of  fluidity  0  and  elasticity  E  in  equation  (8).  The  coefficients 
are  c^  and  Cr. 

(b)  The  EUieoloeical  Parameters  in  the  Creep  Equation 
The  proposed  model  contains  structural  coefficients  and  rheological 
parameters.  Define  creep  compliance  as  creep  per  unit  stress.  Employing 
the  linear  viscoelastic  model  for  creep  (Figure  1),  the  creep  compliance 
for  each  specimen  under  sustained  stress  is  obtained  and  the  rheological 
parameters  evaluated.  Table  III  presents  the  parameters  for  equation  (2) 
and  the  linear  viscoelastic  model.  The  resulting  creep  compliance  (eg) 

--^+   2    1.  (1  -  e  "''"''"')  (9) 


X  00  ,     ^ 

n=l 


n 


The  first  term  in  the  creep  compliance,  containing  the  steady  state 
(ultimate)  viscosity  of  the  concrete  may  be  used  directly  to  obtain  the 
ultimate  fluidity  (0oo  )  required  for  equation  (7). 

The  remaining  Kelvin  elements  form  the  basis  for  determination  of  the 
remaining  parameters  in  equation  (8).  Since  there  will  be  several  (usually 
three  are  sufficient)  Kelvin  elements  in  the  creep  conqiliance,  it  is 
reasonable  to  assume  that  one  of  these  elements  may  be  used  as  the  single 
linear  element  comprising  the  elastic  term  of  equation  (8).  It  may  not  be 
obvious,  however,  which  element  represents  the  delayed  elastic  component. 
Therefore  all  Kelvin  terms  should  be  investigated  by  the  rate  process  theory. 

Initial  fluidity  (0q)  may  be  determined  when  it  is  known  which 
rheological  elements  contribute  to  the  inelastic  portion  of  the  creep. 
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That  is,  after  it  has  been  determined  which  rheological  element  will  be 
used  to  represent  the  delayed  elastic  component,  the  remaining  Kelvin 
terms  of  the  creep  compliance  will  be  used  to  determine  0^^.  The  strain 
rate  function  for  two  or  more  Kelvin  elements  may  be  obtained  from  the 
f  creep  compliance  by  differentiating  with  respect  to  time,  thus  for  Kelvin 

elements  1  and  3 

-0lEit  -0,E^t 

eg  =  0^  e  '^^  ^  +03   e  ^  2  .  (9b) 

At  t  =  0,  ce  o  ~  ^'o'  ^y  definition,  therefore 

00  =  (01  +  03)  .  .  ("> 

C.  The  Analysis  of  Rheological  Assemblies  from  a  Statistical  Mechanics  Viewpoint 

In  the  previous  section  it  was  proposed  that  a  rate  process  study  of 
each  Kelvin  element  in  the  creep  compliance  be  made.  Such  a  study  will 
provide  an  insight  into  the  stress  dependent  nature  of  each  element  considered 
as  a  rheological  assembly.  Each  assembly  contains  the  characteristics  of  a 
molecular  structure  having  energy  levels  representing  the  bonds  between 
similar  assemblies  and  spacings  between  equilibrium  positions  of  the  assemblies 
which  regulate  the  rate  at  which  activation  occurs.  The  rate  process  equations 
presented  in  Chapter  IV  are  therefore  applicable  for  this  rheological  structure. 

For  the  elastic  condition  there  is  expected  a  uniform  statistical  spacing 
between  equilibrium  positions  (Figure  2).  This  assumption  is  based  on  the 
premise  that  for  any  concrete  the  molecular  structure  is  randomly  oriented 
and  elasticity  implies  a  certain  rigidity  of  the  crystalline  structure, 
therefore  the  material  would  be  expected  to  exhibit  uniform  elastic 
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behavior  throughout  the  specimen.  By  the  same  argument  the  free  energy 
of  activation  should  also  be  uniform  throughout  the  specimen.  In  cases 
where  no  specimen  variation  i*  introduced  to  form  comparisons,  such  as 
was  done  in  this  experiment  with  regard  to  shape,  it  may  be  difficult  to 
separate  the  elastically  retarded  element  from  the  other  elements.  In 
some  cases  it  may  be  assumed  that  the  intermediate  activation  energy  is 
related  to  the  elastic  element.  This  may  not  always  be  true,  however, 
where  severe  elastic  deterioration  is  present.  In  general,  however, 
elastic  deterioration  will  require  the  largest  energies  and  viscous  flow 
will  require  the  smallest  energies  of  activation  leaving  intermediate 
levels  of  activation  energy  for  assignment  to  elastic  characteristics. 

The  analysis  of  the  Kelvin  elements  employing  the  statistical 

9 

mechanics  approach  is  derived  generally  from  Herrin's  work  and 

10 
Majidzadeh's  work   with  asphalts.  The  application  of  statistical 

mechanics  to  portland  cement  concretes  is  straight  forward.  The  rheological 

equation  for  the  Kelvin  elements  is  known.  For  each  element  the  strain  is 

f       -0«E„t 
e  -  1;  (1  -  e  '^^  "  )    .  (12) 

It  is  required  to  obtain  an  expression  of  the  element  behavior  in  terms 
of  the  rate  theory.  According  to  the  rate  theory  the  strain  rate,  e  , 
is  a  function  of  stress. 

e  »  A  slnh  B  f  jj  ^  (13) 

where  f-  is  the  decaying  stress  resulting  from  activation  of  rheological 
assemblies  into  new  equilibrium  positions.  For  each  Kelvin  element  we  have 

f  -  f,;  +  f  „  (14) 


\ 
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and 

S     D 
The  total  stress  and  strain  for  an  element  is  f  and  e  respectively.  The 
subscripts  S  and  D  refer  to  the  spring  and  dashpot  components  of  each 
Kelvin  element.  From  equations  (14)  and  (15)  we  may  obtain 

fo  -  f  -  e  E^   .  (16) 

Differentiating  the  strain  with  respect  to  time  for  the  Kelvin  elements 
represented  by  equation  (12)  yields 

-0r,Er.t 

"     .  (17) 

Using  calculated  values  from  equations  (16)  and  (17),  equation  (13)  may 

now  be  solved  by  simultaneous  application  of  least  squares  and  successive 

9 
approximations  for  the  coefficients  A  and  B.  Herrin  has  prepared  a  computer 

program  for  this  solution  which  was  employed  in  this  investigation. 

In  evaltjating  these  data  it  was  considered  that  the  Kelvin  elements 

represented  four  distinct  effects  of  the  creep  behavior  : 

1 .  an  elastic  response 

2..  a  viscous  response 

3.  a  deterioration  or  separation  of  elastic  bonds 

4.  a  structural  growth  of  strength  and  viscosity  due  to  hydration  and  aging. 
Since  these  responses  represent  behavior  of  the  entire  concrete  mass  rather 
than  responses  of  the  separate  phases  comprising  the  concrete,  it  may  be 
assumed  that  the  free  energy  of  activation  and  the  equilibrium  spacing  are 
uniform  for  each  rheological  assembly  (representing  each  Kelvin  element) 
which  is  taken  as  the  basic  assembly  or  flow  unit  in  the  rate  process. 
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Table  IV  contains  the  computed  values  for  the  coefficients  in  the 
rate  function  for  each  type  of  rheological  assembly  in  the  specimen  under 
sustained  stress.  The  method  used  to  determine  these  coefficients  is 
described  and  illustrated  in  Appendix  C.   Some  immediate  relationships 
are  evident,  and  are  described  in  the  following. 

(a)   Spacing  of  Equilibrium  States 

Within  each  concrete  type  the  spacings  between  states  of  equilibrium 
are  relatively  independent  of  shape  for  each  rheological  assembly. 

Furthermore  in  this  investigation,  it  appears  that  these  spacings  are 
also  independent  of  type  of  concrete.   Since  only  the  coarse  aggregate  type 
and  water-cement  (w/c)  ratio  differed  between  these  two  concretes  it  seems 
that  aggregate  and  w/c  ratio  would  have  a  compensating  influence  on  the 
spacing  of  equilibrium  states.  A  small  influence  is  evident  for  rheological 
assembly  number  one,  which  represents  the  long  term  deterioration  effects, 
wherein  the  less  dense  structure  of  the  lightweight  aggregate  yields  a 
slightly  greater  spacing  than  in  the  denser  limestone  concrete.  In  like 
manner,  rheological  assembly  number  three,  representing  short  terra 
viscous  action,  indicates  slightly  greater  spacing  for  the  limestone 
concrete  which  may  be  attributed  to  the  greater  w/c  ratio  in  the  limestone 
concrete. 

The  effects  of  these  two  factors,  w/c  ratio  and  aggregate  strength, 
also  compensated  each  other  from  a  structural  standpoint  to  yield  nearly 
equal  concrete  strengths  which  has  been  considered  by  some  authors  to  be 
an  important  factor  in  creep. '^'  The  major  factor  in  determining  the 
equilibrium  spacing  of  rheological  assemblies  appears  to  be  related  to  the 
concrete  structure  as  measured  physically  by  strength  but  it  appears  that 
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the  most  important  factor  governing  strength  is  related  to  the  thermodynamic 
bonds  between  rheological  assemblies  in  the  mix. 

(b)  Free  Energy  of  Activation 

The  most  significant  influence  on  the  physical  behavior  of  the  concrete 
is  the  level  of  free  energy  of  activation  (height  of  the  potential  energy 
barrier)  for  each  rheological  assembly.   This  level  determines  the  extent 
of  the  thermod3mamic  bonds  and  the  frequency  of  activation  for  the 
contribution  to  creep  from  each  rheological  assembly.   It  can  be  expected 
that  deterioration  effects  resulting  from  permanent  separation  of  elastic 
bonds  will  require  most  energy  and  that  viscous  flow  would  require  least 
energy  for  activation. 
1.  The  second  rheological  assembly. 

The  second  rheological  assembly  representing  elastic  behavior  is  only 
slightly  influenced  by  shape  as  it  affects  hydration  during  the  first  twenty 
days  approximately.  After  about  the  first  month  shape  has  no  effect  on  the 
change  in  the  free  energy  of  activation  in  the  limestone  aggregate  concrete. 
However  shape  appears  to  have  a  slight  influence  on  the  change  in  free 
energy  of  activation  in  the  lightweight  concrete  (refer  to  Table  V).  The 
slightly  larger  average  value  of  the  coefficient  A  from  Table  IV  for  the 
lightweight  concrete  is  indicative  of  a  slightly  greater  elastic  contribution 
than  for  the  normal  concrete.  Therefore,  a  greater  delayed  elastic  recovery 
may  be  expected  from  the  lightweight  concrete  than  from  the  limestone  concrete. 
In  Section  D  elastic  recovery  is  calculated  for  the  specimens  in  test  series  A. 
The  lightweight  specimen  display  a  greater  recovery  than  the  limestone  concrete 
specimen. 
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For  the  second  theological  assembly,  representing  elastic  behavior 
of  the  concrete,  the  results  of  the  statistical  mechanics  approach  shown 
in  Tables  IV  and  V  indicate  the  nature  of  the  stability  of  the  elastic 
structure.  In  the  limestone  concrete  the  first  twenty  days  produce  an 
energy  level  which  is  influenced  by  shape.  However,  following  the  first 
twenty  days  shape  has  no  apparent  influence  on  the  increase  in  energy  as 
seen  from  Table  V.  These  effects  appear  to  be  related  to  the  shrinkage 
characteristics  of  the  concrete  specimen. 

In  the  lightweight  concrete  the  Influence  of  shape  extended  beyond 
the  first  twenty  days.  This  effect  was  probably  caused  by  the  moisture 
contribution  from  the  coarse  aggregate  to  continued  hydration.  In  spite 
of  this  continuing  supply  of  water  in  the  lightweight  aggregate  concrete 
and  a  greater  quantity  of  cement,  the  lightweight  concrete  was  apparently 
deficient  in  moisture  to  achieve  its  potential  activation  energy  level. 
The  limestone  concrete  gained  energy  at  a  greater  rate  than  the  lightweight 
concrete.  The  limestone  concrete  also  contained  a  greater  amount  of  water 
in  the  mix,  though  not  an  excessive  amount,  from  which  a  greater  amount 
of  water  was  colloidally  contained  by  chemical  bond  in  the  gel.  It 
appears,  therefore,  that  the  rate  of  strength  gain  may  be  related  to  the 
quantity  of  free  mix  water  in  the  concrete.  That  is,  even  though  water 
may  be  included  in  the  mix  as  mechanically  contained  moisture  in  the  coarse 
aggregate,  only  a  portion  of  this  water  will  become  chemically  activated 
in  the  hydration  process  after  the  water  is  released  from  the  aggregate. 
The  remainder,  and  perhaps  the  largest  portion,  will  be  treated  as  capillary 
moisture,  ultimately  becoming  associated  with  shrinkage  characteristics  of 
the  concrete  in  moderate  to  small  size  specimens.  In  large  specimens  a 
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greater  portion  of  this  moisture  originating  in  the  coarse  aggregate  may 
become  chemically  hydrated  resulting  in  strength  gain. 

I£  the  level  of  energy  for  the  elastic  state  of  the  structure  is 
considered  to  be  directly  related  to  the  strength,  a  study  of  the 
difference  in  energy  levels  in  the  two  concretes  for  the  circular  shaped 
specimens  would  be  informative.  The  circular  shaped  specimens  were 
identical  in  cross  section  to  the  strength  test  cylinders.  The  limestone 
aggregate  specimen  showed  a  larger  energy  by  about  10  per  cent  from  Table 
IV.  The  limestone  aggregate  specimen  indicated  a  greater  strength  by  about 
12  per  cent  from  strength  tests.  These  results  appear  to  represent 
reasonably  acceptable  trends  but  not  enough  data  are  presented  to  Justify 
any  hypothesis  relating  strength  with  activation  energy  level  of  the 
selected  elastic  rheological  assembly. 
2.  The  third  rheological  assembly. 

The  third  rheological  assembly  represents  short  time  deformation  and 
exhibits  a  decidedly  different  behavior  for  the  normal  limestone  concrete 
than  for  the  lightweight  concrete.  In  the  former,  the  rate  coefficients  A 
are  large  indicating  a  small  activation  energy  resulting  from  the  high 
water-cement  ratio  leading  consequently  to  greater  fluidity.  Uniformity 
of  the  activation  energy  levels  for  each  shape  means  that  the  fluidity 
of  the  system  is  not  influenced  by  the  escape  of  free  capillary  water 
resulting  from  shrinkage  nor  is  the  quantity  of  water  molecules  held  in 
chemical  bond  influenced  by  the  movement  of  free  water  as  it  evacuates 
the  concrete. 

On  the  other  hand,  the  lightweight  concrete  exhibits  a  significant 
influence  of  shape  on  activation  energy.  For  greater  flow  paths  the  energy 
required  for  activation  decreases  implying  greater  fluidity.  Since  it  was 
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seen  that  £ree  water  loss  from  the  cement  gel  does  not  Influence  fluidity 
or  the  free  energy  level,  It  Is  evident  that  moisture  discharged  from  the 
lightweight  aggregate  contributes  to  this  variation.  Evidently,  as 
moisture  Is  discharged  from  the  coarse  lightweight  aggregate  as  free  water 
and  moves  toward  the  surface  of  the  concrete,  portions  of  this  free  water 
are  chemically  absorbed  along  the  way.  With  longer  flow  paths,  greater 
quantities  of  water  are  chemically  absorbed,  causing  greater  fluidity  to 
result  thereby  producing  a  lower  activation  energy.  This  absorption  of 
free  water  would  also  be  expected  to  occur  in  the  normal  concrete, 
however,  the  rate  of  evacuation  of  the  capillaries  is  considerably 
reduced  in  the  lightweight  concrete  due  to  the  slow  emission  of  water  from 
the  aggregaf:e.  This  retarded  moisture  activity  is  evident  in  the  shrinkage 
curves  (Figures  12,  13,  and  14).  The  absorption  of  water  released  by  the  ^ 
aggregate  may  also  be  considered  to  increase  the  water-cement  ratio  for 
the  concrete.  The  significance  of  this  absorption  is  evident  in  Table  V 
from  the  observed  reduction  in  activation  energy  by  about  30  per  cent  in 
the  lightweight  specimen. 

The  third  rheologlcal-assembly  is  changing  during  the  early  life  of 
the  structure.  As  seen  from  Table  V  practically  no  change  occurs  in  the 
limestone  concrete  after  twenty  days.  This  activity  means  that  a  steady 
state  condition  is  achieved  early  in  so  far  as  creep  associated  with  flow 
is  concerned.  Viscous  activity  is  extended  in  the  lightweight  concrete 
due  to  the  reserve  moisture  supply  held  in  the  coarse  aggregate. 

It  appears  reasonable  to  assume  that  specimens  of  larger  sizes  or 
which  are  maintained  in  high  humidity  environments  will  exhibit  this  viscous 
state  through  a  longer  period  of  time  and  possibly  throughout  their  life. 
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Under  these  conditions  the  simple  Kelvin  element  from  which  the  rheological 
assembly  is  developed  may  be  bifurcated  according  to  a  period  of  true 
fluidity  and  a  longer  period  of  quasi  fluidity.  The  former  state  represents 
the  condition  between  initial  pouring  of  the  concrete  and  setting  of  the 
rigid  gel  structure.  The  latter  state  represents  fluid  action  of  the 
elements  still  in  a  colloidal  state  imposed  upon  the  restrictive  nature  of 
the  gel  structure. 
3.  The  first  rheological  assembly. 

The  first  rheological  assembly  exhibiting  long  time  strain  activity 
is  associated  with  large  values  of  activation  energy  which  appear  to  be 
equally  influenced  by  shape  for  both  concretes.  The  more  rapid  evacuation 
of  free  moisture  from  capillaries  of  the  gel  in  the  concrete  specimen  with 
greater  surface  area-to-volume  ratios  results  in  a  reduction  of  the  required 
energy  of  activation  due  to  a  reduction  in  the  rate  of  increase  of  chemical 
bonding  forces.  This  effect  is  directly  related  to  the  amount  of  hydration 
associated  with  each  shape. 

The  larger  values  of  the  coefficient  A  for  the  lightweight  concrete 
represent  a  lower  activation  energy  which  results  from  weaker  bonds  and 
less  hydration.  Weaker  bonds  can  be  attributed  to  the  gel  in  the  lightweight 
concrete  and  to  the  less  dense  Solite  aggregate.  The  gel  may  be  contributing 
to  this  weakness  because  of  insufficient  mix  water  to  hydrate  a  sufficient 
quantity  of  cement  to  form  a  strong  and  rigid  gel  structure.  However,  the 
rate  coefficient  A  also  contains  a  factor  L  (equation  (3c))  related  to  the 
displacement  resulting  from  activation  of  an  assembly,  which  may  be  considered 
to  be  a  modifying  factor  for  non  uniformity  of  behavior.  Stress  concentrations 
and  so-called  plastic  zones  contributing  to  permanent  deformations  result  from 
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piling-up  of  rheological  assemblies  at  internal  obstructions  that  are 
present  within  the  gel  matrix.  Such  obstructions  may  be,  for  example, 
dense  coarse  aggregate  particles  and  reinforcement.  The  conditions  under 
which  deformations  of  this  type  occur  require  suitably  oriented  obstructions 
in  the  path  of  the  migrating  assemblies  with  sufficiently  greater  density 
and  crystalline  rigidity  to  prevent  penetration  of,  or  displacement  by 
the  migrating  assemblies.  This  piling-up  of  assemblies  results  in 
bridging  between  hard  coarse  aggregate  particles  which  transfers  the  stresses 
throughout  the  specimen.  In  soft  aggregate  concrete  the  gel  matrix  is 
more  uniformly  stressed  and  bridging  will  not  be  well  developed  if  it  is 
developed  at  all.  Consequently,  the  strains  resulting  from  activation  in 
the  lightweight  concrete  are  greater  than  in  normal  concrete  where  assemblies 
which  have  created  fully  developed  bridging  prevent  large  deformation  when 
activation  occurs  in  regions  between  the  bridges. 

For  each  of  the  assemblies  discussed  it  must  be  remembered  that 
hydration  will  influence  the  level  of  energy  of  activation  and  perhaps 
the  equilibrium  spacing  also.  Hydration  will  change  the  properties  of 
flow  units  and  adjust  free  energy  levels  causing  formation  of  bonds  thereby; 
affecting  the  potential  energy  barrier  with  time.  Table  IV  contains 
coefficients  A  and  B  in  the  rate  function  for  the  first  twenty  days  only. 
These  values  may  be  compared  with  corresponding  final  values  at  1000  days 
from  Table  IV  for  an  illustration  of  the  influence  of  energy  changes  beyond 
the  first  twenty  days  on  creep  rates.  When  little  change  occurs  in  the  level 
of  activation  energy  for  a  rheological  element  the  steady  state  condition 
is  said  to  exist  for  that  element.  Achievement  of  a  steady  state  condition 
does  not  preclude  further  variations,  however,  if  conditions  in  the  concrete 
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structure  are  changed  so  as  to  cause  a  variation  in  the  quantum 
energy  level.  For  example,  the  lightweight  concrete  experienced  changes 
in  energy  due  to  a  moisture  transfer  from  the  coarse  aggregate  to  the 
gel.   In  this  case  the  increase  in  fluidity  can  be  ascribed  to  the 
discharge  of  free  internal  moisture  from  the  Solite  aggregate.  If 
the  moisture  had  not  been  discharged  into  the  gel,  a  steady  state 
condition  would  have  been  expected  to  occur  sooner. 

Rheological  assemblies  one  and  two  experienced  increases  in 
activation  energy  beyond  twenty  days  as  a  consequence  of  hydration. 
Equilibrium  spacings  in  general  change  proportionately  with  changes  in  the 
free  energy  of  activation  except  in  the  first  rheological  assembly.  The 
large  increase  in  this  case  probably  represents  the  effect  of  greater 
separation  of  molecules  resulting  from  activation  as  the  material  loses 
viscous  damping  capacity  and  elastic  strain  energy  is  converted  directly 
to  kinetic  energy  upon  rupture  of  the  bonds. 

D.  Application  of  the  Structure  Parameters  for  the  Rheological 
Assemblies  to  the  Rheological  Model 

(a)  Retarded  Elastic  Recovery 

Development  of  equation  (8)  representing  the  model  for  the  creep 
mechanism  may  now  be  completed  based  on  the  interpretations  of  the  rate 
process  theory  discussed  in  the  preceding  sections.   Sincfe  the  first 
and  third  rheological  elements  contained  stress  and  time  dependent  vari- 
ations of  non-recoverable  creep  they  are  to  be  included  in  the  second,  non- 
linear term  of  equation  (8).  The  second  rheological  element  satisfies 
elastic  conditions,  and  it  will  therefore  comprise  the  third  term  of 
equation  (8) . 
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A  check  of  the  elastic  Influence  of  this  element  may  be  accomplished  by 
comparing  its  response  with  the  complete  relaxation  data  for  the  actual 
specimen. 

The  strain  accumulated  in  the  elastic  rheological  assemblies  may  be 
computed  at  any  time 


xn  r B  (f •  fp)  ^Asi^iiiA    B^f'-f„)        1 

^^    L      1  '  1  ;     ■"      3- 31        +     5.5:        +.  .  .J. 


This  equation,  derived  from  the  rate  process  theory  for  a  rheological 
element,  is  developed  in  Appendix  A. 

Table  VI,  comparing  results  obtained  from  equation  (18)  by  pooling  all 
specimens  into  lightweight  and  normal  concrete  categories  against  actual 
test  data  supports  the  use  of  the  second  Kelvin  element  to  represent 
retarded  elastic  behavior.  The  calculated  elastic  recoveiry  should  be 
related  to  the  activation  energy  level  at  the  time  of  removal  of  stress. 
For  a  greater  period  under  sustained  loading  the  activation  energy  increases 
thus  the  factor  A  decreases,  thereby  reducing  the  recoverable  elastic  strain. 
The  approximated  final  value  of  the  coefficient  A  is  determined  by  assuming 
that  the  1000  day  value  from  Table  IV  is  the  average  of  the  initial  value 
of  A,  assumed  equal  to  the  20  day  value  from  Table  IV,  and  the  desired 
final  value. 

(b)   Total  Change  in  Fluidity  of  the  Cnncrefg 

Earlier,  initial  and  final  values  of  the  fluidity  terms  0^  and  0oo 
were  determined  from  a  consideration  of  the.  basic  rheological  equations. 
From  a  rate  theory  standpoint  the  fluidity  at  any  time  is 

de 
0  =  df  =  AB  cosh  Bf      .  (19) 
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It  was  shown  in  the.  previous  section  that  the  free  energy  of  activation 
increases  with  time  due  to  the  change  in  the  gel  structure.  Therefore 
the  values  of  A  and  B  in  the  rate  process  equation  are  average  values 
predicting  the  behavior  to  any  arbitrarily  selected  time  Just  as  the 
rheological  parameters  are  average  values  for  any  applied  stress.  If  we 
set  f  =  0  in  equation  (19)  then 

0o  =  AB.  <20) 

which  gives  values  equal  to  the  results  from  equation  (11)  when  the  average 
A  and  B  values  at  any  age  are  used.  However,  any  value  of  f  greater  than 
zero  will  result  in  greater  fluidity  than  0^  when  employing  equation  (19). 
Therefore  it  becomes  obvious  that  the  fluidity  change  must  include  the 
true  variation  in  internal  energy  and  equilibrium  spacing  for  an  accurate 
account  of  the  non  recoverable  strains  at  any  stage  of  the  creep  process. 

An  understanding  of  the  nature  of  the  change  in  the  quantum  states  of 
the  concrete  is  not  within  the  scope  of  this  study.  The  true  value  of  0^ 
should  be  approximately  equal  to  the  value  determined  from  equation  (20). 
For  the  statistical  mechanics  approach  which  is  represented  by  equation  (20), 
the  values  of  0^  =  AB  do  not  change  when  using  the  data  at  twenty  days  or 
at  1000  days  from  Table  IV. 

Equation  (11)  provides  a  decrease  in  fluidity  with  respect  to  time 
whereas  equation  (20)  provides  an  increase  in  fluidity  with  respect  to  stress. 
The  true  fluidity  is  probably  a  non  linear  function  of  time  and  stress 
requiring  a  knowledge  of  the  physico-chemical  thermodynamic  variations 
affecting  the  quantum  states  of  the  mix. 

Assuming  the  ultimate  fluidity  0oo  to  be  established  by  equation  (10)  the 
total  change  in  fluidity  may  be  approximated  by  assuming  a  variation  between 
the  values  established  by  equations  (11  or  20)  for  0^  and  equation  (10)  for  0bo 
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The  total  change  In  fluidity  is  influenced  primarily  by  the  mix 
proportions  of  the  concrete  and  bj-  the  ambient  conditions  of  the 
environment  in  which  the  concrete  will  be  used.  For  most  structural 
concretes  in  any  environment  the  gel  structure  becomes  sufficiently 
rigid  after  a  long  time  that  little  or  no  viscous  creep  occurs  even  in 
the  presence  of  high  humidity.  Therefore  it  may  be  justified  to  use 
zero  as  a  value  for  0oo  ,  If  this  expedient  step  is  taken  it  only  becomes 
necessary  to  evaluate  0^  from  test  data  obtained  from  the  concrete  under 
initial  loading. 

(c)  Inelastic  Behavior  and  the  Cnefficlenta  nf  SttrurMiral  R^flhn■^^^v 
The  coefficients  of  structural  stability  for  equation  (7)  may  be 

determined  by  evaluating  the  responses  from  the  control  specimens, 

2    2 
-(ct  +  Cf  (fo  -  f  )) 
0  =  000  +  (0o  -  0«.)  e  .  (7) 

From  a  consideration  of  the  behavior  of  the  control  specimens  (Series  A) 
for  which  f  -  f^,  a  determination  of  the  time  dependent  component  in 
the  structural  stability  coefficient  may  be  made  by  solving  for  the  aging 
component,  Cj..  Substituting  f  =  f q  into  equation  (7)  yields  the  solution 
for  the  time  dependent  component  of  creep  only, 

0  -  (iloo "^t  (21) 


0^.0 


e 


Taking  the  log  of  both  sides  results  in  a  solution  for  c^. 

0  -  0«, 
c  =  -  In  :    "  (22) 


in  which 

€vis 


(23) 
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^  vis  ^^   simply  obtained  by  subtracting  shrinkage,  Initial  elastic  strains 
and  delayed  elastic  strains  from  the  total  creep  strain  for  each  specimen. 
Time  dependent  coefficients  (c.)  for  the  non-recoverable  creep  component 
are  presented  In  Figures  10  and  11.  The  coefficient  c^  Is  assumed  to  be 
a  function  of  time  only,  therefore  It  Is  entirely  Independent  of  stress 
variations.  Equation  (22)  Is  used  to  obtain  the  value  of  c^^  for  the 
doncrete  from  experimental  data.  However  equation  (22)  may  not  be  used 
as  a  general  relation  for  the  prediction  of  creep  behavior.  The  coefficient 
c^  must  Include,  In  addition  to  the  variable  time,  all  of  the  Influencing 
factors  which  contribute  to  the  aging  characteristics  of  creep  of  concrete. 
Hansen  and  Mattock   have  evaluated  shape  effect  on  creep  and  presented 
a  relationship  between  creep  and  shape  In  an  exponential  form.  It. Is  not 
Incongruous  to  assign  a  similar  relationship  to  other  factors  which 
Influence  creep.  Accordingly  the  coefficient  c^  is  assigned  an  exponential 
form 

ct  =  t"  e    "     ,  (24) 

where  F^(t)  is  a  polynomial  function  of  time.   The  equations  c^  are  given 
in  Table  VII  for  each  control  specimen.  The  constants  in  the  c  function 
are  related  to  the  time  dependent  factors  influencing  concrete  aging, 
hydration  for  example,  which  cause  variation  of  the  energy  states  in  the 
concrete  Independently  of  stress.  The  polynomial  exponent  has  been  limited 
to  a  cubic  function.  The  limitation  is  based  partly  on  an  acceptable  fit 
of  the  function  to  the  experimental  data  as  represented  by  the  error  sum 
of  squares  in  Table  VII.  The  closeness  of  fit  is  indicated  in  Figures  10 
and  11.  The  plots  were  made  from  a  computerized  curve  fitting  program. 
The  exponent  was  limited  to  a  cubic  function  also  because  it  was  quite 
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obvious  that  higher  power  terms  would  only  tend  to  influence  the  creep 
at  advanced  ages  and  in  not  a  highly  significant  manner.  A  reasonably 

successful  fit  may  be  obtained  in  the  early  ages  by  substitution  of  a 

16 
constant  for  the  polynomial  function  as  was  done  by  Hansen  and  Mattock. 

The  complete  form  of  the  aging  component  c^.  of  the  structural 

stability  coefficient  is 

c,  -t^e-<^l  +  ^2'  +  ^3t^"'^^'^^   .  (24a) 

From  Table  VII  the  constant  a  is  relatively  independent  of  shape  or 
concrete  mix.  Further,  it  is  reasonably  close  to  unity  so  that  setting 
8=1  does  not  result  in  a  serious  error  for  the  concrete  specimens 
tested.  The  polynomial  exponent  is,  however,  related  to  shape.  By 
employing  the  circular  shape  as  a  standard,  since  it  has  a  weighted  flow 
path  length  equal  to  unity  for  a  six  inch  diameter  specimen,  all  other 
shapes  or  sizes  may  be  related  to  the  six  inch  diameter  cylinder  as  a 
standard.  Redesignating  the  coefficients  b^^  in  the  polynomial  function 
to  Cjj  for  the  six  inch  diameter  cylinder,  the  coefficients  b^  may  be 
approximately  determined  from 

FK  (1  -  w)    "I 

n=l 

where  w  is  the  weighted  flow  distance  for  the  concrete  specimen  from 
Table  II  and  K  is  a  constant  associated  with  the  concrete  mix.  K  is 
equal  to  six  for  the  lightweight  Solite  concrete  used  in  this  experiment 
and  K  is  equal  to  three  for  the  limestone  concrete. 
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The  stress  dependent  component,  c,,  in  the  structural  stability 
coefficient  is  obtained  from  the  control  specimens  employing  the  rate 
coefficients  from  the  rate  process  theory  presented  in  Table  IV.   In 
Appendix  A,  the  creep  for  a  single  rheological  assembly  was  determined. 
The  combined  strain  for  two  such  assemblies  representing  the  Kelvin 
elements  one  and  three  from  the  rheological  model  of  Figure  1  may  be 
written  in  the  form  of  an  infinite  series, 

e,  +e  o  =  2   [  (ABi  X  /E)i  +  (AB^  \  /E),]  Ilfl_:_l^> 

j=l,2^  ^  •  ^  •      (26) 

where  i  =  2j  -  1. 

The  component  C£  is  related  to  stress  changes  in  the  concrete 

beyond  the  initial  applied  stress  f  .   If  the  stress  remains  unchanged 

then  only  aging  influences  the  creep  rate  and  the  C£  component  is  not 

effective.  However  under  any  change  in  stress,  increase  or  decrease,  the 

creep  rate  will  be  influenced  according  to  the  effect  of  the  stress  on  the 

energy  required  for  activation  of  a  rheological  assembly.   If  a  creep  test 

is  performed  on  a  standard  cylindrical  specimen  and  the  macroscopic  creep 

data  evaluated  by  viscoelastic  methods,  the  resulting  rheological  solution 

for  the  viscous  non  Newtonian  elements,  such  as  Kelvin  elements  one  and 

three  in  this  investigation,  may  be  converted  to  the  form  represented 

by  equation  (26) .   This  process  results  in  an  evaluation  of  non  Newtonian 

concrete  strain  in  terms  of  stress  change  and  the  linear  viscoelastic 

parameters  from  the  rheological  study.  The  accuracy  of  the  strains 

determined  from  equation  (26)  is  limited  to  the  accuracy  of  the  rheological 

parameters  0^^  and  E^j  in  reflecting  the  true  concrete  behavior.   This 
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limitation  on  accuracy  is,  however,  directly  related  to  the  degree  of 
work  which  the  analyst  is  willing  to  expend  in  obtaining  more  accurate 

viscoelastic  parameters.   Surely  i£  methods  as  described  by  Freudenthal 

4 
and  Roll  are  employed,  the  greater  effort  will  yield  better  results. 

The  important  aspect  of  the  relationship  presented  in  equation  (26) 
is  in  having  a  function  which  is  sensitive  to  changes  in  stress  rather 
than  time.   This  relationship  makes  it  possible  for  the  structural  analyst 
to  evaluate  the  influence  of  stress  changes  on  creep  behavior  with  greater 
precision  than  was  possible  from  just  a  rheological  study. 

In  order  to  relate  these  non  Newtonian  strains  to  the  component  c^ 
of  the  structural  stability  coefficient,  let 


N. 


(AB^  \   /E)^  +   (AB^  \  /E)3 


Equating  the  strains  represented  by  equation  (26)  from  the  average  states 
of  the  rheological  assemblies  one  and  three  to  the  stress  dependent 
creep  term  in  equation  (8)  for  the  rheological  model  yields 

(CJ  -0oo)  exp  I  -Cf  (f/  -  f^)   1  ft  =  2     Ni  <y  '  V       (27) 


Dividing  by  stress  f  and  time  t,  substituting  from  Appendix  A 

t  =  -  ^  In  I- 
E     ^o 


(A-4) 


and  taking  the  log  of  both  sides  of  equation  (27),  the  stress  component  of 
the  structural  stability  coefficient  is  determined 


2    2) 
Cf  (f„  -  f'^=ln 


r     r             1  (f'"' 

I             (AB^)!  +  (AB^)3|           i 
1=1.2 

-  ^9  /^>    1 

(0o-0eo)      In  f/fo 

(28) 
»i=2j-l. 
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The  coefficients  A,  B,  0^  and  0oo  have  also  been  determined.  Therefore 
knowing  the  stress  variation,  the  right  side  of  equation  (28)  may  be  substituted 

^f  (^0^ 
from  equation  (8), 


2    2 
for  the  coefficient  cr  (f   -  f  )  and  the  total  creep  may  be  calculated 


(8) 
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CHAPTER  VII 


ANALYSIS  OF  RESULTS 


A.  Application  of  the  Model  to  Concrete  Under  Sustained  Stress  (Test  Series  A) 

The  model  is  designed  from  the  response  of  test  series  A  in  which  the 
concrete  specimen  were  maintained  under  sustained  stress  for  over  1000  days. 
Hence  it  is  expected  that  the  model  will  represent  the  behavior  well. 
Figure  18  confirms  the  comparison  between  creep  data  and  model  for  each 
condition  of  specimen  shape  and  concrete  mix. 

The  creep  data  for  Figure  18  were  obtained  by  subtracting  the 
corresponding  shrinkage  from  Figures  12,  13,  or  14  from  the  total  strain 
for  the  specimen  shown  in  Figures  15,  16,  or  17,  with  allowance  made  for 
shrinkage  occurring  up  to  the  time  of  stressing.  There  is  a  prominent 
seasonal  effect  on  shrinkage  and  total  strain  which  becomes  inconspicuous 
in  the  carry-over  to  the  creep  data.  The  extent  of  the  seasonal  effects 
in  the  creep  data  when  compared  with  the  effects  In  shrinkage  is  a  measure 
of  the  interaction  between- creep  and  shrinkage  of  concrete. 

The  model  which  has  been  developed  is  in  the  general  form  represented 
by  equation  (8) . 

f   r       ^  -(Ct.  +  Cf(fo  -  f^)K     f  r     -0eEe«=l 

The  elastic  components  are  not  considered  to  be  Influenced  by  shape. 
Therefore  E,  0e>  Eg  are  the  rheological  parameters  determined  directly 
from  the  vlscoelastlc  analysis  of  the  data  from  the  control  specimen. 
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The  inelastic  component  is  a  function  of  shape.  The  four  parameters 
involved  are  0oo  ,  0q,  Cj.  and  c..  Each  parameter  appears  to  be  influenced 
by  shape  or  flow  path.  The  value  for  0eo  may  be  obtained  from  equation  (10) 
employing  the  results  in  Table  III. 

«„  -il     ;^-"<«=-"^  (29) 

where  0^     is  the  ultimate  fluidity  for  the  standard  control  specimen,  the 
six  inch  diameter  cylinder,  and  w  is  the  weighted  flow  distance  from  Table  II 
for  the  specimen  under  consideration.  Equation  (29)  is  applicable  to  either 
concrete.  When  the  standard  six  inch  diameter  cylinder  is  used,  W^  =  1.0 
and  equation  (29)  becomes 

The  value  for  0o  is  obtained  from  equations  (11)  or  (20).  For  the 

specimens  tested,  employing  values  from  Table  III  in  equation  (11)  results 

in  no  noticeable  influence  of  shape  on  0^  for  the  limestone  concrete. 

However,  the  lightweight  Solite  concrete  exhibited  a  variation  in  0^ 

influenced  by  shape.  The  reason  for  this  difference  in  behavior  may  be 

attributed  to  the  proximity  of  the  moisture  laden  Solite  aggregate  to  the 

concrete  surface.  The  apparent  difference  in  viscosity  is  caused  by 

differences  in  moisture  loss  from  the  aggregate.  Employing  equation  (11) 

from  page  41, 

-1.45  M  (1  -  w)  ,_. 

00  -  (0i  +  8>3)c  e  .  (30) 

where  the  subscript  c  refers  to  the  values  for  a  standard  control  specimen, 
a  six  inch  diameter  cylinder  in  this  case. 
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M  =  0  for  standard  limestone  concrete 

M  =  1  for  Solite  lightweight  concrete  as  used  in  this  experiment. 
A  determination  of  the  variation  of  M  is  not  within  the  scope  of  this 
experiment . 

The  structural  stability  coefficients  c.  and  c^  have  been  fully 
described  in  the  previous  chapter.   The  aging  component  c^.  was  defined 
by  equations  (24a)  and  (25).   For  the  concrete  in  this  investigation  the 
value,  of  c.  was  obt.iined  by  omitting  the  power  of  t  in  equation  (24a). 

,^  =  t  e-^'l  "  '2^^  ■"   '3''  ■"  '^''>         .  (31) 


where  h^  =  C^ 


K  (1  -  w) 
^    7+  1 

Ln=l   ^n 


(25) 


The  values  of  C^  are  determined  by  fitting  equation  (31),  with  C^   replacing 
b_,  to  the  creep  data  from  a  standard  control  specimen  made  from  the 
concrete  under  investigation.  The  values  for  c  are  determined  for  any 
time  t  by  use  of  equation  (22).   At  least  ten  time  intervals  should  be  used 
in  obtaining  t.hfi  values  of  C^^.   Values  for  b^  for  any  other  shape  or  size 
of  specimen  are  obtained  from  equation  (25)  where  K  =  3  (1  +  M) 
K  =  3  for  limestone  aggregate  concrete 
K  =  6  for  Solite  aggregate  concrete. 

The  value  for  the  structural  stability  component  Cf  is  determined  by 
substitution  defined  by  equation  (28).   Values  for  0^  and  0«)   are  obtained 
as  described  earlier.  Only  the  rate  parameter  A  is  dependent  upon  flow  path 
or  shape  (see  Table  IV).   Determination  of  the  values  for  the  parameter  A 
to  be  used  in  equation  (28)  may  be  approximated  from  values  obtained  for  a 
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standard  control  specimen  of  six  inch  diameter  for  each  inelastic  rheologlcal 
assembly  unit.  For  the  first  rheologlcal  assembly 

^  ,  ^^,  ^(1  -  K/lOXl  -  w)     _  (3„ 

where  A,^  is  the  calculated  value  of  A  for  the  first  rheologlcal  assembly 
of  the  cylindrical  control  specimen.  K  and  w  have  been  previously  defined. 
The  value  for  the  third  rheologlcal  assembly  is 

A3  =  A3=  .-'■'  «<'-">       .  (33) 

where  M  has  been  previously  defined. 

All  parameters  may  be  obtained  from  a  single  sustained  creep  test   on 
a  standard  control  specimen.  However,  in  each  of  the  relations  above  it  is 
not  possible  to  cross  from  one  concrete  mix  to  another.  Therefore  the  test 
must  be  made  on  the  type  mix  to  be  used  in  the  prototype  structure.  The 
limited  scope  of  this  experiment  precludes  development  of  a  model  which  will 
include  mix  variations  and  environmental  variations  in  its  prediction  of 
creep  behavior. 

An  analysis  of  the  model  (Figure  9)  reveals  certain  characteristics 
of  its  behavior.  Considering  the  first  and  third  terms  of  equation  (8)  it 
is  evident  that  these  components  of  the  model  reflect  no  variations  in  the 
molecular  structure.  The  variations  due  to  Irrecoverable  viscous  flow, 
destruction  of  bonds,  and  growth  of  new  bonds  are  therefore  all  contained 
in  the  second  term  of  equation  (8).  In  the  case  of  concrete  under  sustained 
stress  these  factors  are  all  time  dependent  and  the  magnitudes  of  their 
influence  on  inelastic  creep  is  directly  proportional  to  the  sustained 
stress  level. 
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Froportlonallty  of  creep  and  stress  Is  limited  to  low  stress  ranges. 
When  the  applied  stress  exceeds  about  30  per  cent  o£  the  twenty-eight  day 

strength,  f,  the  creep  rate  Increases  and  proportionality  between  creep 

4 
and  stress  ceases.  Freudenthal  and  Roll  have  determined  that  only  the 

Inelastic  components  of  creep  are  not  proportional  to  stress.  In  the 

model  represented  by  equation  (8)  therefore  only  the  second  term  representing 

Inelastic  creep  components  would  be  non  linear  with  respect  to  stress  for 

high  stress  levels.  This  Investigation  did  not  Include  a  study  of  stress 

4 
level.  Freudenthal  and  Roll  Investigated  stress  level  on  different  mixes 

producing  concretes  varying  from  very  viscous  to  very  fluid  In  nature. 
Employing  their  results  equation  (8)  may  be  modified  to  Include  consid- 
eration of  non  proportional  creep  response  under  high  stress  conditions. 

Since  the  second  term  of  equation  (8)  Is  linear  with  respect  to  stress 
level  any  creep  component  obtained  from  It  at  any  time  would  be  the  linear 
creep  component  only.  In  order  to  obtain  the  true  creep  the  linear  creep 
component  must  be  multiplied  by  a  factor  which  produces  the  non  linear 
component  of  creep.  The  non  linear  component  Is  related  to  the  viscosity 
of  the  mix.  A  more  viscous  mix  will  yield  a  smaller  non  linear  creep 
component.  Since  Freudenthal  and  Roll's  results  represented  the  extreme 
conditions  of  viscosity  It  Is  not  known  how  the  non  linear  creep  component 
varies  with  viscosity.  Therefore  It  Is  assumed  that  the  non  linear  creep 
component  Is  Inversely  proportional  to  the  Initial  viscosity  of  the  mix  at 
the  time  the  stress  Is  applied. 

In  corporatlng  non  proportional  relations  between  creep  and  stress  Into 
the  model,  equation  (8)  becomes 
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C  =  |f  [0»-.  (0,-0<„)  e-^^t+'^f  C^o'-^'))]  (1  +p)  ft  +  1^  [l  -e-«»«'«^J.(8a) 
where  p  =  (2.5  +  0^  x  10^)  (f/f^  -  .30) 

The  ratio  f/f  represents  the  ratio  of  stress  level  to  strength  for  concrete 
being  evaluated.  It  is  tacitly  assumed  here  that  the  stress  to  strength 
ratio  at  which  proportionality  ceases  is  about  0.30.  The  value  for  p  may 
not  be  negative.  Therefore,  the  value  for  p   is  set  equal  to  zero  for  all 
stress  levels  resulting  in  values  for  f/f^  less  than  0.30. 

While  grouping  all  inelastic  effects  in  the  second  term  is  mathematically 
convenient  there  are  certain  disadvantages  with  respect  to  the  activity  of 
the  third  term.  When  evaluating  the  delayed  elastic  recovery  upon  unloading 
the  concrete  the  elastic  state,  or  more  precisely  the  quantum  state  of  the 
molecular  structure  at  the  time  of  unloading  determines  the  amount  of 
recovery.  The  constant  coefficients  of  the  third  term  are  unable  to 
satisfactorily  predict  this  variation.  This  deficiency  may  also  affect  •.. 
recovery  values  for  concrete  under  decreasing  stress.  Further  investigation 
is  required  to  fully  understand  the  manner  in  which  the  elastic  response 
coefficients  vary. 

B.  Application  of  the  Model  to  Concrete  Under  Decreasing  Stress  (Test  Series 
B  and  C) 

The  model's  coefficients  were  derived  from  the  response  of  the  specimen 
held  under  sustained  stress  for  convenience.  The  application  of  the  model 
to  the  specimens  of  test  series  B  (stress  decrease  of  65  per  cent)  and  C 
(stress  decrease  of  50  per  cent)  is  therefore  independent  of  the  test  results 
for  these  specimens.  Figure  19  compares  the  model  with  creep  test  data  for 
specimens  of  test  series  B.  The  model  agrees  with  the  data. 
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Figure  20  compares  the  model  with  creep  test  data  for  specimens  of 
test  series  C.  The  agreement  Is  not  as  good  as  for  test  series  B.  The 
apparent  cause  of  the  dlscrepency  seems  to  be  related  to  the  delayed 
elastic  element  (third  term)  which  recovers  a  greater  amount  of  strain  than 
It  should.  This  effect  was  described  earlier  In  connection  with  delayed 
elastic  recovery  of  the  specimen  of  series  A.  Not  only  was  the  average 
value  of  the  structure  coefficient  A  used  In  connection  with  elastic 
recovery,  but  the  coefficient  A  was  also  declared  Independent  of  time  in 
the  derivation  of  the  stress  dependent  factor  for  the  stability  coefficient. 
These  two  assumptions,  based  on  an  Insufficient  knowledge  of  the 
thermodynamic  properties  of  the  material  with  time,  lead  to  low  predicted 
values  of  concrete  strain.  If  It  were  known  how  the  coefficient  A 
decreases  with  time  and  mix  then  a  smaller  recovery  could  be  predicted 
and  accuracy  Improved.  This  effect  only  influences  creep  predictions  at 
advanced  ages.  The  model  appears  to  adequately  predict  creep  in  the  early 
ages  under  stress.  The  early  influence  of  fluidity  on  creep  is  controlled 
to  a  great  extent  by  the  third  Kelvin  element  of  the  elementary  model  of 
Figure  1.  This  element  Is-not  greatly  influenced  by  time  beyond  20  days 
(Table  V).  Therefore  the  assumption  of  constancy  of  thermodynamic 
coefficients  A  for  this  rheologlcal  assembly  does  not  introduce  serious  error. 

The  response  of  the  model  for  test  series  B  agrees  better  than  the 
response  for  test  series  C  because  the  error  caused  by  the  dlscrepency  in 
the  third  term  of  equation  (8)  is  proportional  to  stress.  Since  series  C 
had  a  larger  average  stress  the  deviation  of  the  model  from  the  data  will 
be  greater. 
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Equations  (6)  or  (8a)  are  treated  like  any  vlscoelastic  equation 
with  respect  to  varying  stress.  Super-position  of  responses  from  each 
term  is  made  for  time  intervals  over  which  the  applied  stress  is  assumed 
to  be  constant.  The  solutions  for  test  series  B  and  C  were  obtained  by 
super-position  of  responses  over  time  intervals  for  which  the  stress 
changed  by  2.0%  of  the  original  applied  stress. 

C.   Influence  of  Shape  on  Creep 

A  study  of  the  structure  constants  in  Table  IV  shows  the  influence 
of  shape  on  the  transient  creep  components.  For  the  size  specimen  used 
in  this  study  shape  has  only  a  slight  influence  on  creep  behavior  in 
limestone  aggregate  concrete.  This  effect  is  associated  almost  entirely 
with  the  action  of  the  first  rheological  assembly.  However,  shape  does 
significantly  influence  the  early  creep  rate  in  lightweight  aggregate 
concrete  of  the  type  used  due  to  the  retarded  shrinkage  effect.  This 
would  lead  one  to  believe  that  a  similar  effect  would  result  due  to  size 
difference.  That  la,  for  sufficiently  longer  flow  paths,  even  in  limestone 
concrete,  where  a  sufficient  amount  of  free  moisture  still  remains  in  the 
concrete  after  initial  hydration  is  completed,  creep  may  be  Influenced 
by  fluidity  variations  similar  to  the  effect  caused  by  shape  in  the 
lightweight  concrete  in  this  experiment.  It  seems  quite  likely  that 
shape  effects  reported  by  some  authors  may  really  be  an  interaction 
between  size  and  shape  in  the  specimen  tested. 

Even  though  early  viscous  creep  in  the  lightweight  concrete  was 
Influenced  by  shape  there  appeared  to  be  no  continuation  of  this  behavior 
into  longer  periods  of  time.  There  was  no  effective  ultimate  creep 
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dlf ference  among  shapes  from  which  it  is  concluded  that  the  long  term 
deterioration  creep  component  was  Inversely  influenced  by  shape  and 
therefore  tended  to  compensate  for  the  initially  higher  creep  in  those 
specimen  with  high  volume -to -surface  area  ratios. 

It  may  be  concluded  that  shape  has  a  binary  effect  on  creep  behavior: 

1.  An  increase  in  bonds  due  to  greater  hydration  in  connection  with  longer 
flow  paths  results  in  reduced  viscous  creep. 

2.  A  reduction  in  fluidity  (resulting  in  greater  bonds)  leads  to  reduced 
viscous  creep.  Fluidity  is  reduced  as  a  result  of  rapid  moisture  loss 
associated  with  shorter  flow  paths. 

Both  effects  are  influenced  by  size  also,  perhaps  to  a  greater  extent. 
These  creep  effects  are  further  altered  by  stress  to  such  an  extent  that 
shape  may  appear  not  to  have  any  Influence  at  all.  Evidence  of  this  behavior 
is  Indicated  in  Table  VIII  wherein  specimens  under  decreasing  stress  did 
not  exhibit  the  same  shape  dependency  as  specimens  under  sustained  stress 
over  a  long  term  period. 

D.  Effect  of  Aggrogate  on  Creep 

In  addition  to  the  interaction  effects  of  aggregate  and  shape  described 
in  the  previous  section  there  was  a  significant  Influence  noted  relative  to 
the  stress  level  after  a  period  of  time.  The  Initial  elastic  modulus  which 
is  dependent  upon  the  density  of  the  concrete  was  naturally  smaller  for  the 
lightweight  concrete  thereby  producing  greater  elastic  strains.  However, 
in  addition  to  this  expected  behavior,  the  lightweight  aggregate  appeared 
to  Influence  transient  creep  in  such  a  way  that  creep  in  the  lightweight 
concrete  was  greater  under  high  stress  and  less  under  low  stress  than  creep 
in  the  limestone  concrete  under  similar  stress  conditions. 
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TABLE  VIII 


TOTAL  OBSERVED  TRANSIENT  AND  STEADY  STATE  CREEP  OF 
CONCRETE  UNDER  VARYING  STRESS  LEVELS 

After  1000  days 


Type  of 
Concrete 


Creep  Creep  Creep 

Shape      Test  Series  A   Test  Series  B   Test  Series  C 
Micro  in/ in     Micro  in/ in     Micro  in/ in 


Limestone 

Circle 

490 

200 

310 

Rectangle 

610 

200 

290 

Concrete 

Cross 

705 

245 

300 

Avg.  3  Shapes 

600 

215 

300 

Lightweight 

Circle 

575 

160 

275 

(Solite) 

Rectangle 

640 

175 

310 

Concrete 

Cross 

730 

180 

280 

Avg.  3  Shapes 

650 

170 

290 
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Two  different  factors  caused  this  effect.   In  a  low  stress  situation 
the  elastic  properties  dominate  the  creep  behavior  of  the  concrete.   It 
was  shown  in  Chapter  VI,  Section  C  that  the  gel  in  the  lightweight  concrete 
was  more  elastic  than  the  gel  in  the  limestone  concrete.  Therefore,  elastic 
recovery  was  greater  in  the  lightweight  concrete  resulting  in  a  smaller 
ultimate  creep.  Under  high  stresses  the  irrecoverable  term  became  dominant. 

Due  to  its  greater  rate  of  irrecoverable  strain,  manifested  by  the 
relative  magnitudes  of  the  coefficients  A  (Table  IV)  for  the  first 
rheological  assembly,  the  lightweight  concrete  experienced  greater  creep 
over  a  long  period  of  time.  This  creep  is  undoubtedly  a  result  of  the 
permanent  surface  distortions  of  the  coarse  aggregate  particles  under 
stress. 

Reference  to  the  comparisons  made  in  Table  VIII  between  average 
transient  creep  for  the  two  concrete  types  provides  a  clear  picture  of  the 
stress  level  effects  on  creep.  Under  high  stress,  test  series  A,  creep 
in  the  Solite  concrete  was  greater  whereas  under  low  stress,  test  series 
B,  creep  in  the  Solite  concrete  was  less  than  in  the  limestone  concrete. 


CHAPTER  VIII 


CONCLUSIONS 


A.     General 

It  has  been  shown  that  creep  behavior  may  be  adequately  predicted  from 
macroscopic  observations  by  employing  structure  parameters  in  a  rheological 
model.  These  parameters  include  the  change  in  fluidity  due  to  natural 
aging  effects  under  stress  and  the  change  in  fluidity  due  to  stress 
variations.  A  rheological  model  may  be  constructed  with  structure 
coefficients  from  statistical  mechanics  theory  which  will  predict, 
reasonably  well,  the  creep  of  concrete  from  a  consideration  of  the  response 
from  a  test  of  several  specimens  sustained  at  low  constant  stress.  The 
parameters  which  are  used  in  this  model  represent: 

1.  The  rate  of  transition  of  the  gel  from  fluid  to  solid. 

2.  The  total  expected  change  in  the  structure  of  the  gel. 

3.  The  purely  elastic  action  of  the  gel  and  aggregate  system. 

The  many  aspects  of  the  influence  of  mix  on  the  viscoelastic  nature 
of  concrete  precludes  development  of  a  generalized  creep  model  at  this  time. 
It  is  proposed  in  this  report  that  an  analysis  of  creep  behavior  for  a 
particular  concrete  is  possible  and  a  prediction  made  from  test  specimens 
of  standard  dimensions.  The  tests  should  be  conducted  on  concrete  from 
the  same  mix  as  the  prototype  structure  and  in  the  same  general  environment. 
Using  equation  (8a)  the  creep  may  be  predicted  at  any  time  under  varying 
stress  conditions  and  for  any  conventional  size  or  shape  of  structure.  Mass 
concrete  structures  do  not  satisfy  conditions  for  use  of  this  model. 
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B.   Test  Results 

Maintaining  equal  paste  contents  and  nearly  equal  early  strength 
characteristics  for  two  concrete  mixes  produced  creep  responses  as  follows: 

1.  Shape  did  not  appear  to  Influence  the  early  creep  of  specimen 
of  the  size  tested  when  rapid  shrinkage  occurred.  If  for  some  reason, 
such  as  moisture  retaining  aggregate  or  greater  volume  for  example,  a 
certain  amount  of  free  moisture  Is  maintained  In  the  gel  beyond  the 
Initial  phase  of  early  hydration  shape  may  become  Influential  as  It 
affects  the  rate  of  change  of  fluidity.  This  creep  response.  In  the  early 
stages  of  hydration,  to  shape  Is  due  largely  to  a  rapid  decrease  In  fluidity 
In  high  early  strength  concretes.  Since  It  Is  a  rate  comparison  of 
fluidity  change  which  Is  being  made  therefore,  normal  strength-gain 
concretes  should  show  a  greater  Influence  of  shape  on  creep  than  high 
early-strength  concretes. 

2.  Coarse  aggregate  has  a  diverse  Influence  on  the  elastic  properties 
of  concrete.  It  Is  already  known  that  density  Influences  the  Initial 
elastic  modulus  (ACI  code,  1963).  An  Increase  In  density  results  In  concrete 
with  greater  elastic  modulus.   However,  density  also  Is  a  measure  of  the 
rigidity  or  stiffness  of  the  aggregate,  and  aggregates  which  are  stlffer 
than  the  gel  structure  appear  to  produce  a  less  elastic  nature  over  a 
long  period  than  aggregates  which  are  less  stiff  for  equal  strength 
concretes.  This  result  stems  from  the  fact  that  the  gel  In  concrete 
with  stlffer  coarse  aggregate  must  have  poorer  quality  than  the  gel  In 
equal  strength  concrete  made  from  less  stiff  coarse  aggregate, 

Whereas  the  gel  behavior  Is  more  elastic  In  the  lightweight  concrete, 
the  long  term  non  recoverable  creep  Is  greater  because  of  more  general 
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permanent  damage  to  the  rough  and  cellular  surface  of  the  softer  aggregate 
by  the  stiffer  gel.   In  concrete  with  harder  aggregate  crushing  of  the 
gel  is  localized  and  less  irrecoverable  creep  results. 

3.  Water-cement  ratio  influences  the  initial  fluidity  and  the 
strength  gain.  Lower  water-cement  ratios  materially  reduce  the  initial 
creep  of  concrete  by  producing  a  lower  initial  fluidity  and  they  reduce 
the  over-all  creep  activity  because  of  an  increased  level  of  activation 
energy. 

4.  The  amount  of  moisture  retained  by  the  coarse  aggregate  in 
lightweight  concretes  has  a  significant  binary  affect  on  the  fluidity 
change  after  initial  hydration.  On  one  hand  the  supply  of  moisture 
during  aging  produces  greater  levels  of  activation  energy  from  increased 
hydration  thereby  reducing  creep  and  on  the  other  hand  the  fluidity  is 
increased  as  moisture  is  absorbed  by  the  gel  matrix  thereby  increasing 
creep.  While  these  two  factors  tend  to  offset  each  other,  over  the  long 
period  either  effect  may  become  dominant  depending  on  the  test  conditions. 
Therefore  the  creep  response  may  take  different  aspects  dependent  upon 
the  quantity  of  moisture  included  in  the  coarse  aggregate  at  the  time  of 
mixing . 

C.  The  Model  for  Creep  Prediction 

In  order  to  take  into  account  variations  of  coarse  aggregate,  the  model 
was  developed  with  a  parameter  M.  M  is  considered  to  represent  a  property 
of  the  coarse  aggregate  of  the  concrete.  Evaluation  of  M,  however,  was 
outside  the  scope  of  this  investigation.   In  order  to  satisfy  the  model  M 
was  found  to  be  equal  to  0.0  for  the  limestone  concrete  used  in  this  study, 
and  it  was  equal  to  1.0  for  the  lightweight  concrete  used.  It  is  presumed 
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that  these  aggregates  represent,  to  ail  intent  and  purpose,  the  limits  of 
coarse  aggregate  structure.  That  is,  for  soft,  porous  aggregate,  such  as 
Solite,  M  achieves  a  maximum  value  of  1,0  and  for  hard,  dense  aggregate, 
such  as  limestone,  M  has  a  minimum  value  of  0.0.  Thus  it  is  presumed  that 
M  varies  between  the  limits  of  0.0  and  1.0  for  all  types  of  coarse  aggregate. 
However  there  is  no  experimental  evidence  to  substantiate  this  assumption. 
For  this  reason  it  is  suggested  that  these  values  for  M  be  applied  only 
to  structural  concretes  with  similar  aggregates  to  those  used  in  this 
investigation  until  verification  can  be  established. 

On  this  basis  equation  (8)  may  be  employed  to  predict  creep  in  concrete 
from  a  creep  test  on  standard  six-inch  diameter  specimens  under  low,  - 
sustained  loading.  The  parameters  in  equation  (8)  may  be  determined  as 
follows: 

000  -  Equation  (29a) 
0Q  -  Equation  (30) 

c   -  Equations  (22),  (25),  (31)  in  consecutive  order 
Cc     -  Equation  (28)  employing  equations  (32)  and  (33). 
The  rheological  parameters  required  to  satisfy  the  above  relations  are 
obtained  by  conventional  rheological  theory  as  described  in  reference  22. 

In  order  to  expand  the  accuracy  of  equation  (8),  results  from  research 

4 
reported  by  Freudenthal  and  Roll  were  incorporated  to  account  for  the  non- 
linear relationship  between  creep  and  stress  at  high  stress  levels. 
Accordingly,  equation  (8a)  may  be  used  when  stress  levels  exceed  0.30  f^. 
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D.  Procedure  for  Using  the  Model 

In  order  to  facilitate  the  use  of  the  model  a  brief  procedural  outline 
is  presented. 

If  it  is  desired  to  understand  the  behavior  of  a  structural  concrete  in 
service  with  regard  to  its  creep  character,  a  test  series  may  be  undertaken 
on  about  four  6'-lnch  concrete  cylinders  under  low  sustained  stress  of 
approximately  20  per  cent  of  f\ 

1.  In  normal  environments  subjected  to  seasonal  variations  the  concrete 
in  the  structure  will  eventually  achieve  a  state  where  the  ultimate  fluidity 

(000)  is  very  small  and  maybe  assumed  equal  to  zero.  Therefore  the  test 
need  only  be  run  on  the  concrete  cylinders  long  enough  to  establish  the  slope 

and  intercept  of  Ross'  parametric  creep  curve.  Perhaps  up  to  90  days  of 

testing  may  be  required.  The  averaged  data  from- the  four  specimens  will  be 

plotted  in  the  form  of  a  straight  line 

t  "  b  +  a  t/e  , 

where  t  is  the  ordinate  and  t/e   is  the- abscissa.  The  intercept  b  and  the 
slope  a  are  obtained  directly,  and  theymay  be  substituted  Into  the  creep 
equation  (1) 

e  -  at/(b  -ft).  (1) 

2.  Having  established  the  basic  parametric  relationship  between  creep 
strain  and  time  a  linear  rheologlcal  model  may  be  developed  as  shown  in 
Figure  1  which  is  represented  by  equation  (2).  The  method  for  determining 
the  parameters  in  equation  (2)  and  the  procedure  for  developing  the  linear 
vlsco-elastic  model  is  presented  in  reference  22,  or  any  treatise  on 
conventional  rheologlcal  theory. 
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3.  The  solution  £rom  step  2  provides  knowledge  of  the  parameters  0  „       and 
E  for  the  final  model  shown  In  Figure  9  represented  by  equations  ^  or  (Ba ) . 
E  is  the  elastic  modulus.  0oo  is  calculated  using  equation  10.  If  the 
preliminary  test  is  not  carried  sufficiently  far  to  evaluate  0oo  .it  may  be 
set  equal  to  zero  with  very  little  probable  error. 

4.  A  rate  process  solution  must  then- be  performed  on  each  of  the  Kelvin 
elements  from  step  2.  Equation  (12)  represents  the  behavior  of  each  Kelvin 
element.  To  represent  the  rate  of  creep  as  a  function  of  stress  in  the  form 
of  equation  (13)  an  iteration  procedure  is  performed.  In  order  to  solve  for 
the  coefficients  A  and  B  in  equation  (13)^  e   and  f^  must  be  determined  for 

each  element  in  accordance  with  equations  (17) and  <t6) respectively  at  corresponding 
times.  Equation  (l3) , 

G  =  A  sinh  B  f  (13) 

may  now  be  solved  by  the  iteration  procedure  outlined  in  Appendix  C.  A 
value  for  the  coefficient  B  is  assumed  and  -substituted  into  each  of  the  two 
equations  for  A  given  in  Appendix  C.  All  values  are  known  for  the  summations 
in  each  equation.  The  summations  are  made  over  time  increments  and  the  value 
of  A  is  calculated  from  each  equation  and  compared.  Each  step  in  the  iters ti^p 
process  is  performed  with  a  diffetent  value  for  B.  When  the  difference  between 
the  calculated  values  of  A  begins  to  deviate  from  zero,  as  shown  in  Appendix  C, 
these  values  of  A  and  B  correspond  to  solution  of  the  rate  equatior  (13). 
The  summations  taken  over  time  increments  are  summed  to  the  time  at 
which  prediction  is  required.  This  process  is  quite  tedious  and  lends  itself 
quite  conveniently  to  computer  solution. 
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5.  The  rheological  assembly  corresponding  to  elastic  behavior  is  selected 
as  described  earlier  in  the  report.  The  remaining  Kelvin  elements  may  then 
be  used  to  determine  0^  from  equation  0.1) 

6.  To  determine  the  stability  coefficient  c^,  the  viscous  strains  (  e^j^g) 
are  determined  from  the  simple  rheological  model  of  figure  1  by  eliminating 
the  elastic  components. 

e  ,  -  ft/v«,    +  2       f/E  (1  -  e"^nV)   . 
vis      '^*       viscous 

The  summation  is  made  over  all  remaining  Kelvin  elements  after  the  delayed 
elastic  term  has- been  removed.  The  coefficient  c^.  for  the  test  cylinder  at 
any  time  is  determined  from  equations  (22!>  and  (23) .  This  procedure  results  in 
a  set  of  values  of  c^.  corresponding  to  Increasing  values  of  time.  Therefore 
a  mathematical  solution  for  c  is  required  in  the  form  presented  by  equation 
(24a)with  c_  termsi  corresponding  to  the  control  cylindei;,  substituted  for  the 

b_  terms. 

7.  In  projecting  the  test  results  to  predict  creep  in  a  prototype  member 
the  shape  and  size' influence  must  be  considered.  The  weighted  flow  path 
distance  (w)  is  a  measure  of  each  effect.  The  weighted  distance  for  the 
6-inch  cylinder  is  equal  to  1.0.  For  any  other  member  weighted  distance 
may  be  calculated  by  dividing  the  member  cross-section  into  small  increments 
of  area,  and  multiplying  these  small  areas  by  their  respective  shortest 
distance  to  the  member's  surface  then  dividing  the  product  by  the  total 

cross-sectional  area.  , 

8.  The  creep  behavior  of  the  concrete  is  influenced  greatly  by  the  coarse 
aggregate.  Only  two  types  of  aggregate  have  been  tested  for  their  influence 
on  creep,  Brooksville  limestone  and  Florida  Solite.  Since  Florida  Sollte 
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is  a  porous  aggregate,  its  capacity  to  bear  moisture  is  great.  Values  for  the 
aggregate  coefficient  M,  representing  the  influence  of  the  aggregate  on  creep 
were  found  to  be  zero  for  the  limestone  aggregate  and  1.0  for  Solite  aggregate 
with  24  per  cent  initial  moisture  by  weight;  It  may  be  presumed  that  the 
value  of  1.0  will  vary  for  different  moisture  contents,  however  no 
experimental  evidence  is  available. 

9.  The  aging  coefficient  of  stability  may  now  be  projected  for  the  member 
response  from'  the  test  results.  •  Employing  equations (£ 5 )and (24a*.,  Cj.  may  be 
evaluated  for  the  prototype  member. 

In  similar  fashion,  0tc   and  0^  may  be  projected  to  the  prototype 
member  as  outlined  in  section  C  of  this  chapter.  Also  rheological  structure 
coefficients  A  for  the  long  term  and  short  term*  transient  conditions  of 
creep  may  be  projected  by  using  equations  62) and  030  respectively. 

10.  The  creep  at  any  time  may  now  be  obtained  for  any  stress  history  by 
employing  the  rheological  model  of  Figure  9  represented  by  equations  (8) or  (3a^ 
Since  the  equation  is  linear  with  respect  to  stress  (f),  in  order  to  obtain 
solutions  under  varying  stress  conditions  piecewise  time  increments  having 
constant  stress  must  be  employed  and  superposition  of  the  effects  must  be 
made.  This  procedure  is  outlined  in  detail  in  reference  22. 

In  employing  equation (8),  the  value -of  the  stress  dependent  stability 
coefficient  (c,)  is  obtained  simultaneously  by  equation  (28) 

If  the  stress  at  anytime  exceeds  0.30  f'  the  non- linear  components  of 
creep  must  be  considered.  Then  equation  (8a)is  applicable  without  variation 
in  the  procedure. 
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APPENDIX  A 

CALCULATION  OF  STRAIN  FROM  THE 
ACTIVATION  OF  A  RHE0L06ICAL  ASSEMBLY 

The  strain  for  one  Kelvin  element  Is 

f        -0nEnt 
e  =  i;  (1  -  e      )    .  (A-1) 

The  strain  rate  written  In  terms  of  the  structure  coefficients  is 

^  =  e   =  A  sinh  B  f.      .  (A-2) 

dt  ^ 

The  stress  fp  may  be  written  in  terms  of  time  by  substituting  equation  (A-1) 

into  equation  (16) 

-0  E  t 
fo  =  f  -  e  E„  =  fe  "  "     .  (A-3) 

Solving  equation  (A-3)  for  time  t  yields, 

J: ,    ^D  ,.  ,v 

Also  differentiating  (A-3)  with  respect  to  time  yields 

j£  -0  E  t 

££d  =  -  f  0^E^e  "  "  (A-5) 

from  which 

1    1     0nEnt 

Rewriting  equation  (A-2)  and  integrating  yields  the  strain 

e  =  A  J*   sinh  B  fp  dt     .  (A-7) 
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Substituting  equation  (A-4)  and  (A-6)   into  (A-7)  produces 

^-  -At    HTv-    J*   f"     9i^h  B  fn  df_  (A-8) 

f     ^n%         fj,  D       D 

which  when  integrated  yields 

BfD    B^fD^   B^fjj^         -T]  .  ■ 

-  (In  fjj  -  TTi;  +  JTl   -  37^7  +  .  •  .)J     +  C    (A-9) 

simplifying  to 

e  =  -  A  TTT-   r  Hfi-  +  lis.  +  ^V  +  .  .  .   "I   +  c  .  (A-lO) 
Initially,  at  t  «  0,  fp  =  f^,  £  =  0,  therefore, 

which  when  substituted  into  (A-10)  yields 

e  =A-i-  Tli^o-i-M  +!ii!2!:J2!  4. B^fo^-  fp^^  .  .   .1  (A-12) 

0nEn  Ll-1!       3'31       S'SI         J 


or  , 

e=Aj^   "    "'  ^'°'  "  'P  ^   >  (A»13) 

''nEn   Z      i  •  i  • 
j=l,2 

where  i  =  2j  -  1. 


APPENDIX  B 

DEVELOPMENT  OF  THE  STRUCTURAL  FACTORS 
FOR  COEFFICIENT  OF  STABILITY 

The  Influence  of  aging  and  normal  stress,  f,  is  known  to  change 
the  flow  characteristics  of  concrete.  The  shear  stress  on  any  element  is 
equal  to  f/2  for  a  uniaxial  stress  f .  According  to  Reiner  ^    fluidity 
is  related  to  stress  squared  (f  /4).   The  coefficient  of  stability  of  the 
concrete  structure  may  be  written  as  a  function  of  the  time,  and  stress 
rate  of  change  of  fluidity  of  the  material.  If  all  influences  on 
structural  stability  are  included  in  this  coefficient  a  general  relation 
will  result.  Fluidity  changes  differently  with  time  under  vairying 
stress  conditions.  Therefore  the  structural  stability  coefficient  must 
be  an  interaction  function  of  tine  and  stress.  Let 

be  the  generalized  form  of  the  structural  stability  coefficient. 

According  to  Reiner  ^   aging  of  the  mix  alters  its  physical 
properties.  Since  the  values  of  all  constants  of  physical  state  of  the 
mix  are  affected,  when  the  numerical  value  of  a  constant  is  given  the 
corresponding  physical  state  must  also  be  given  with  regard  to  temperature 
and  pressure.  If  these  are  not  mentioned,  one  may  assume  that  the  constant 
is  changing  with  the  physical  states  of  the  mix.  Therefore  X  may  be 
assumed  to  vary  with  time  and  we  may  therefore  separate  the  structural 
stability  coefficient  into  time  and  time-stress  components. 
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If  the  structural  stability  coefficient  X  is  composed  of  two  factors, 
one  a  time  dependent  factor  c^  related  to  aging,  and  the  other  a  factor  c^ 
related  to  stress  dependent  variations  with  time,  then  the  general  form  of 
the  structural  stability  coefficient  may  take  the  form 

i=-^  ^  +-a!£f  .  ^  S^   ,  (B.2) 

where  c^  is  the  aging  component  of  the  structural  stability  coefficient, 
dependent  on  time  only. 

Cf   is  the  stress  dependent  component  of  structural  stability, 
dependent  on  time  and  stress. 

The  negative  signs  are  applied  because  X  increases  with  increasing  time 
whereas  X  decreases  with  increasing  stress. 

The  boundary  conditions  for  equation  (B-1)  are: 

at  t  =  0  ,    0  =  00  and  f  =  fg  (B-3) 

att=",    ^"S'eo  approximately  .  (B-A) 

For  the  boundary  conditions  (B-3)  it  is  assumed  that  loading  is  accomplished 
before  the  viscosity  changes  by  a  significant  amount.  In  the  application 
of  the  methods  for  determining  viscosity  of  the  mix,  a  load  must  be  applied 
to  the  concrete.  Therefore  the  initial  viscosity  is  assumed  to  be  the 
viscosity  at  the  time  load  is  applied. 

The  solution  of  equation  (B-1)  is  achieved  by  rewritting  in  the  form 

5  ^  (0  -  00°  )  =  .  S(f^/4)St  .B-5^ 

(0  -  000  )  X  .  •  ■    ^ 

The  solution  of  equation  (B-5)  is 

0  =  0^  +  Ae"^*=  e"°f  ,  (B-6) 
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where  Cj^  is  a  function  of  time  only  and  Cr  is  a  function  of  stress  and 

time.  Applying  boundary  conditions  (B-3)  results  in  a  solution  for  A 

A  -  00  -  0eo      .  (B-7) 

This  solution  tacitly  implies  that  the  components  Cf.  and  Cf   are 

proportional  to  time  and  stress  as  follows: 

Ct  =  c^  (t)  (B-8) 

Cf  =  cf  (f^^  -  f^)          .  (B-9) 

Boundary  condition  (B-4)  is  satisfied  for  any  stress  change  condition. 

The  equation  for  the  change  in  fluidity  under  varying  time  and  stress 

conditions  is  written 

-  (  Ct(t)  +  Cf  (f^  -  fS) 
0  =  000  +  (00  -  000  )  e  .  (B-10) 

Increases  in  stress  over  the  initial  stress  is  seen  to  increase 
fluidity  and  creep,  whereas  decreases  in  stress  from  the  initial  applied 
stress  decreases  fluidity  and  creep.  Since  linear  proportionality  is 
implied  by  the  component  c^^  the  parameter  t  is  removed  from  equation  (B-10), 
Therefore 

0  =  000  +(00-  0bo  ^,\(-f^C,if\-fh)       ^  ^^_^^^ 


APPENDIX  C 

METHOD  FOR  DETERMINATION  OF  THE  COEFFICIENTS 
IN  THE  RATE  PROCESS  EQUATION 

e  =  A  sinh  Bf 
The  computer  program  developed  by  Herrin  and  used  by  Majldzadeh 
determines  the  value  of  the  coefficient  A  from  the  known  initial  stress 
and  complete  strain  data  for  a  Kelvin  element.  The  value  of  A  calculated 
from  each  of  the  two  solutions 

S  '^  e  "  "  sinh  Bfo 

A,  =  -J^ (C-1) 

2 
2   sinh  Bfjj 

and 

„  f       -(J„E„t 
t  r  e       fo  cosh  Bf D 

A2  = (C-2) 

2  J.  fjj  sinh  Bfj)  cosh  Bfp 

are  compared  for  successively  approximated  values  of  B.  Solution  for  the 
values  of  A  and  B  is  declared  at  the  iteration  step  in  which  the  compared 
values  of  A  begin  to  steadily  deviate  from  one  another.  The  development 
of  equations  (C-1)  and  (C-2)were  by  M,  Herrin,  Department  of  Civil  Engineering, 
University  of  Illinois,  Urbana,  Illinois. 

Figure  C-1  illustrates  the  typical  mode  of  deviation  of  the  A 
values  for  each  iteration  step  and  the  method  of  solution  of  the 

transient  creep  associated  with  each  rheological  assembly  for  a  single 

specimen. 
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